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Abstract 


This study examines the evolution of the Klein-Nishina (KN) cross section for hydrogen 
and platinum atoms in relativistic and non-relativistic scenarios. Using online MATLAB 
tools, the constructed model of KN cross section for both situations is calculated. The 
observation demonstrates that as electron energy increases, electronic and atomic KN 
cross section lowers, which lessens the relationship between the incoming electron and 
the target atom. The importance of atomic size is demonstrated by a comparison of the 
KN cross sections of hydrogen and platinum, with platinum possessing a higher cross 
section because of its greater size and more scattering centers. Higher atomic numbers 
result in greater interaction between sites and bigger cross sections, and they also play a 
role in this relationship. The variations in cross sections are also influenced by electron 
binding effects, especially in the instance of hydrogen. It is demonstrated that relativistic 
effects have a major impact on the KN cross section, including an exponential decline 
seen as the energy of electrons increases. In the relativistic domain, the electron’s 
increased mass and decreased wavelength increase scattering probability and produce a 
greater cross section than in the non-relativistic case. This research sheds light on the 
variables affecting KN cross sections and their influence on radiation therapy, medical 
imaging, and materials science. 
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1. Introduction 


In recent years, fuel cells made from hydrogen have drawn interest as a 
practical alternative energy source. Generally speaking, fuel cells use fuels 
like the gas hydrogen, methanol, or natural gas to transform the chemical 
energy contained in the fuel they contain into direct current (DC) without 
intermediary the combustion process, as is the case of traditional engines [1]. 
Fuel cells are a potential candidate for an environmentally friendly form of 
energy technology in the future due to their increased efficiency and low 
emissions. The fuel cells market and, consequently, the shift to a hydrogen 
economy are determined by the kind of fuel cell, the cost of the components, 
and the production processes. The price of fuel cells has dropped by 80% in 
the last couple of decades (since the 1990s), thanks to technological 
developments. Additionally, it is anticipated that the price of platinum will 
decline due to the rising demand for natural resources [2]. An fuel cell is an 
electrochemical machine that uses two redox processes to wansform the 
chemical energy of a fuel (typically hydrogen) plus an oxidizing agent 
(commonly oxygen) into electrical energy [3]. Over the past three decades, 
as clean energy technologies are sought after globally for environmental and 
sustainable development goals, the proton exchange membrane fuel cell 
(PEMFC) has attracted a lot of interest. PEMFCs are devices which employ 
hydrogen fuel and transform chemical energies directly into electricity, 
bypassing the limitations of the Carnot cycle. As a result, they are more 
efficient and emit almost no emissions in comparison to conventional 
electricity-generating technologies and direct consuming of fossil fuels. 


A single PEMFC is typically built in a hierarchical form, with a proton 
exchange membrane serving as the center layer through which primarily 
protons can pass. The catalyst layer (CL) of the anode and cathode, the gas 
diffusion layer (GDL) between the two electrodes, and the bipolar plates are 
atranged in order from the center to the edges. Together, the membrane, CLs, 
and GDLs are referred to as membrane electrode assemblies (MEA), and 
their physical makeup could be compared to that of porous media with multi- 
scale pores. Grooves called gas flow channels (GFCs) are mechanically 
carved into the innermost layer of the polar plate [4]. During cell operation, 
hydrogen and air are delivered into the cell via the flow channel and, in turn, 
pass through the anode's and cathode's respective gas diffusion layers to reach 
the catalyst layer. The molecules of hydrogen are depleted of their electrons 
and broken down into protons in the catalyst, which is the layer of the anode. 
To create an electric current, electrons are required to move through an 
external circuit. The formed electron, proton and residue goes interaction in 
effective area called, scattering cross sections. The scattering region can be 


ee 


calculated by using equation (1) from differential cross section phenomenon 


[5]. 


do v2 (k'\" 
Eras (;-) (1 + cos*6) (1) 


In terms of the wavelengths of the photons, Compton described this 
relationship as follows: 


A’ -2A=1-cosd (2) 


where A and A’ are the photon's pre- and post-deflection Compton 
wavelengths [6]. 


In this work authors are applying Klein-Nishina (KN) cross section in the 
framework of electrode scattering which helpful insights about the 
connection among electrons and the electrode materials. In order to maximize 
the effectiveness and effectiveness of PEMFCs, it is essential to comprehend 
the KN cross section in this particular system. Electrochemical reactions in 
PEMFCs are greatly aided by electrode materials like platinum or other 
catalysts. As they engage with electrode surfaces during fuel cell operation, 
produced electrons experience scattering events. The likelihood of each of 
these scattering events occurring is measured using the KN cross section. The 
information regarding scattering in PEMFCs and its effects on performance 
can be examined by the KN cross section around an electrode scattering. On 
taking account the KN cross section in electrodes dispersion within PEMFCs 
offers useful knowledge for enhancing the design, functionality, and 
longevity of these energy-converting devices. It may direct the choice of 
materials and electrode engineering, and it can promote PEMFC technology 
for a variety of uses, such as transport and stationary power generation. 


II. Literature Review 


The initial gas-based voltaic cell powered by a mixture of hydrogen and 
oxygen was created by Sir William Grove in 1839, the same year he also put 
forth the notion of the fuel cell's effect. Ostwald laid the theoretical 
groundwork for fuel cell research in 1894. When attempting to understand 
the rate of electromechanical processes in fuel cells, Van Santen developed 
the Ostwald-proposed step rule concept [7]. In the 1920s, Lawaczek created 
numerous schemes for hydrogen-powered automobiles, trucks, buses, trains, 
and engines. He also made a substantial contribution to the development of 
pressurized electrolyzes. Although knowing that the idea for fuel cells first 
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emerged in the beginning of the 1840s, it wasn't until Francis T. Bacon began 
researching them between 1932 and 1952 that this technology was once again 
developed, allowing the very first time the development and evaluation of a 
5 kW hydrogen fuel cell stack. Later, this technology was created for the life 
support, guiding, and navigation systems used in NASA's Gemini and Apollo 
programs [8]. During the course of the last 20 years, the fuel cell sector has 
made enormous strides, leading to discussions about making investments in 
hydrogen. The development of fuel cell automobiles was a focus for several 
decades for automakers from various nations. The first commercially 
available fuel cell vehicle was released by Toyota in 2014, capping years of 
research and development work. As a result, the public began to view fuel 
cell cars as a significant component of the next generation of mobility rather 
than:'as an experimental technology. Countries including China, the US, 
Japan, and a number of European nations concentrated their efforts on 
advancing this technology throughout the course of the following five years 


[9]. 


From 2022 and 2027, the worldwide fuel cell industry is anticipated to 
increase from an USD 3 billion economy to a USD 9 billion economy at a 
26.0% CAGR [10]. Over the course of the last few years, scientific research 
has also grown more interested in fuel cells. The peak growth in patent 
applications in the hydrogen-based fuel cell industry [11] (23% increase) 
from 2016 to 2020 confirms the aforementioned claim. The development of 
these technologies is driven by personal goals and international obligations 
in nations including the United States of America, China, India, Japan, South 
Korea, and several European nations. On the other hand, the efficiency of the 
fuel cell is affected by the quantum species (electron, proton, hydrogen and 
oxygen) formed during the interaction of inlet fuel (hydrogen and oxygen) 
with platinum. KN cross section during the collision of free electrons with 
atoms (H-atom and Pt-atom) close to the cathode of proton exchange 
membrane fuel cells (PEMFCs) is the goal of this work. During the collision 
of free electrons with H and Pt atoms, respectively, the maximum KN cross 
section measured for single scattering is approximately 70.2 m* and 66.0 m2. 
The maximum KN cross section recorded for 1ml flow of hydrogen is about 
—26.6 m? and —22.3 m? [12]. In 2004, Roy and Pratt used X-ray generator to 
test KN cross section theory with observations (from 11 to 40 keV) of the 
entire atom and the observations shows conventional sources ranged from 5% 
to 50%, but the measurements made with a synchrotron source indicated a 
small percentage difference. The Compton scattering measures in a range 
below 10 keV with a broader range of works using angles other than 90 and 
distinct atomic number Z values in order to confirm the sufficiency of the 
theoretical treatment. 
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The KN cross section behavior for H and Pt atoms is the main topic of the 
current investigation. Future studies may look into a larger range of topics to 
develop a more thorough grasp of the elements affecting the KN cross 
section. Understanding how these variables affect the KN cross section would 
require research into many elements with various atomic sizes, atomic 
numbers, and material characteristics. Dhobi et al. conducted research on the 
KN cross section for non-relativistic scenarios and the scattering zone 
surrounding the electrode in 2023. Furthermore, Dhobi et al. (2022) did not 
analyze the KN cross section in either the relativistic or non-relativistic case, 
but rather the interaction of quantum species created around the electrode of 
PEMFC [13]. There are numerous projects underway to increase the 
PEMFCs cell's efficiency through experimentation, theory, modeling, etc. 
Additionally, study is- being done on scattering fields theoretically, 
practically, virtually, etc. However, there aren't many PEMFC and scattering 
coupled works. Given that PEMFC and scattered are both novel in the fields 
of research, academia, and industry. Consequently, it is crucial to research 
dispersion around the PEMFC electrode. 


III. Methods and Materials 


Hydrogen fuel cell electric vehicles (FCEVs), when compared with 
traditional internal combustion engine vehicles, they are more efficient from 
a different perspective. They have a range for driving of more than 300 miles 
and can be fueled similarly to vehicles with traditional internal combustion 
engines in around 5 minutes. The technology of FCEVs boost efficiency, like 
regenerative brakes that collect the electrical power lost during brakes and 
store it in batteries [14]. Due to their ease of use, quiet operation, high 
efficiency, and modular design, fuel cells (FCs) have attracted a lot of 
attention in the automotive sector in recent decades. The usage of FCs in 
electric automobiles (EVs) is predicted to grow quickly, spark a revolution, 
and eventually replace conventional vehicles, according to technological 
breakthroughs. Commercial cars, programs, and research demonstrate that 
efforts are being made to guarantee that FCEVs have performance 
improvements sufficient for daily transportation requirements. First, the 
application areas, distinguishing characteristics, and operational 
circumstances of FC types and electrical motors are reviewed. According to 
their structure, complexity, and architectural frequency of use, power 
converters—voltage regulating and motor drive topologies—that are 
employed in FCEVs are described in depth [15]. 


-6- 
Free electrons, protons, and leftover hydrogen are created when hydrogen 
comes into touch with the electrodes' catalyst. When the formation electron 
collides with the atom of hydrogen, a great deal of collisions occurs in the 
system, which causes heat to be produced. However, in this instance, we take 
into account the collision that occurs of the hydrogen atom and a free 
electron. PEMFCs have a serious issue with fuel cell efficiency because of 
the heat produced during collisions. We are interested in studying the Kelin- 
Gordan differential cross section that results from the collision of an electron 
and a hydrogen atom. The classical Thomson scattering can be described 
using a generalized Klein-Nishina differential equation for Compton 
scattering with a limited train of pulses. The quantity of scattering photons 
produced at a particular solid angle dQ is 


ANscat _ ie Egc|Ex(w)|? do ‘ (3) 
dQ 2nh|w| dO 


= CO 


And because the scattered photon has energy hw’, the total scattered energy 
is: 


dU, “ egclE,(w)|? w' dao 

<= | oC lEx(w)| sd (4) 
dQ, is 21 w dQ, 

The input gamma ray impacts an atomic electron, causing atomic ionization, 
and this is the primary scattering mechanism that produces 511 Kev photons. 
The KN differential cross section formula will determine the angle at which 


the incident photon will scatter, 


de\ _. Zrg 1 2 2 a*(1-—cosd)? 
(3), me (seat coass) (« eos" 8) steer) (5) 


A variation of the differential atomic cross sectional area equations for K-N 
can be found here. Additionally, Knoll discovered the complete KN cross 
section per atom in 1989; 


=e, ((@) mode 6 
Og = 20 A do) (6) 


Where @ is scattering angle overall photons. Now from (5) and (6), we get: 


my es 
WZTe 1 : 
= Et (eer eee 6 
Oa an | + (scar) (Cr ) 


a?(1 —cos6)? sate ; 
- [1 + a(1 — cos6)] on ) 


After solving, the following is the total KN cross section per atom: 


2 2 fie pote, move In(i+2a@) _ _1+3a 
Pa = Zante a2 | 1+2a a + 2a ao} (8) 


Since electronic cross-sections are multiplied by each element's charge 
number Z to get Klein-Nishina atomic cross-sections, equation (7) states that 
0, = Z.¢, the electronic cross-sectional area for KN is 


= 2 {tte pare _ me In(1+2@) —_1+3ar 
Oe = 2nTe a2 1 1+2a a = 2a aaa} (9) 
where r, = 2.818 fm is the classical electron radius, Z is the nuclear charge 


= — = ———,, for non relativistic case. For 
of the target molecule and a mec? OSLiMeV” 


eae _ me? mec? 
relativistic case, energy of electron E, = Taye ie 


0.511MeV [17]. Also, E, = Te therefore for relativistic case a, = 


Now equation (11) can be modified as 


po 2 ieee Pater) marten In(l+2ar) ae 
Oe = 21Te az | 14+2a, ay - 2ar (1+2a,)2 (10) 


Also from, dg = Z.a, we have 
re 2 (itary 2(1+ar) a In(1+2a@,) In(i+2a,) 3: ae} 
Og = 22ZnTe a2 | 1+2ay Oy |+ 2a (1+2a,)2 (11) 


TV. Results and Discussion 


[16] and mc? = 


1 
J1—v2 ce? 


The KN cross sections for electron and atom interactions are shown in Figure 
1 as a function of electron energy. With increased electron energy, the cross 
section is seen to decrease. The cross section is initially large at lower energy, 
indicating a greater likelihood of interaction. But as electron energy rises, the 
cross section gradually contracts until it reaches a fixed value. It is clear from 
comparing the cross sections of platinum and hydrogen that platinum has a 
higher KN cross section than hydrogen and electronic interaction. The size of 
the atoms under consideration accounts for the majority of this variation. 


R's 


Platinum atoms’ increased size increases the scattering likelihood and raises 
the cross section. 


By taking into account the energy dependency of scattering processes, it is 
possible to explain the decreasing pattern of the KN cross section with rising 
electron energy. Because their energy scales are more closely matched at 
lower energies, the relationship between an incoming electron and the object 
being targeted atom has a greater impact. The contact, however, weakens with 
increasing electron energy, resulting in a decrease in the width of the cross 
section. The significance of atomic size is demonstrated by comparing the 
cross sections of hydrogen and platinum. The cross section of platinum is 
higher because platinum atoms are bigger than hydrogen atoms. This is 
explained by the fact that there are more scattering centers accessible, which 
increases the likelihood of scattering occurrences. The KN cross section 
obtains a constant value in the high-energy area. This tendency implies that 
the incident electron's interaction with the target atom weakens and becomes 
less reliant on electron energy. Due to the predominance of relativistic effects, 
whereby the wavelength of the input electron shrinks dramatically relative to 
the size of the scattering centers, the interaction is diminished. 
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Figure 1: KN cross section with electron energy 


Figures 2, 3, and 4 depict the KN electron and atomic cross section together 
with the electron's velocity. The KN electronic and atomic cross sections fall 
exponentially with increasing energy of electrons in the relativistic situation. 
The two main causes of this phenomena are the contraction of the electron 
wavelength and the relativistic mass increase. For the relativistic situation of 
mass, an electron's mass grows noticeably when it approaches relativistic 
energies, which are energy close to the speed of light. The relativistic mass 


nO 
of an object traveling at a high speed is determined by, in accordance with 
Einstein's theory of special relativity. The trajectory of the electrons during 
scattering events is impacted by the increase in mass. At greater energies, the 
electron behaves resembles a heavy particle because its relativistic mass is 
much more than its rest mass. The KN cross section decreases as a result of 
the higher mass since it lowers the likelihood of scattering. 


In the relativistic wavelength situation, the electron's wavelength gets shorter 
as its energy rises. An electron's de Broglie wavelength is calculated using 
the formula (A = h/p), where A is the wavelength, h is Planck's constant, and 
p is the electron's momentum. A shortened de Broglie wavelength results 
from the electron's momentum increasing as its energy does. By decreasing 
the true dimension of the scattered centers (atoms), the diminishing electron 
wavelength has an impact on the scattering process. The efficiency of 
scattering decreases when the wavelength of the electron approaches or falls 
below the dimensions of the scattering centers. The KN cross section 
decreases as a result of this decline in scattering effectiveness. As a result, the 
KN electronic and atomic cross sections exponentially drop with increasing 
electron energy due to the interaction between the relativistic mass increase 
and the decreasing electron wavelength in the relativistic regime. 
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Figure 2: KN Electronic cross section with electron energy 
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Figure 3: H-Atom KN cross section with electron energy 
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Figure 4: Pt-Atom KN cross section with electron energy 


The KN cross section for Pt is greater than the KN cross section for hydrogen 
(H), as well as the KN electronic cross section, when the KN cross sections 
from Figures 2, 3, and 4 are compared. Additionally, it is discovered that the 
KN cross section for hydrogen is greater than that of the KN electronic cross 
section. There are various underlying causes for this variation in cross section 
values. The atomic size of the particles that interact is a significant factor for 
the observed variations in the KN cross sections. Hydrogen atoms are much 
smaller than platinum atoms. Due to their bigger size, platinum atoms offer a 
larger target region for the incoming electrons to engage with during the 
scattering process. Since there are more scattering events as a result, platinum 
has a greater cross section than hydrogen and electronic interactions. 


a2 


The target atoms' atomic number is another element that affects the cross 
section discrepancy. Compared to hydrogen, Pt has a higher atomic number. 
Higher atomic numbers equate to more electrons in the atom's outer shells, 
which increases the number of potential places for interaction for the 
incoming electrons and affects the scattering likelihood. As a result, platinum 
has a greater cross section than electrical and hydrogen interactions. The 
difference in cross section between electronic interactions and hydrogen is 
also influenced by the electron binding effects. The attractive force of the 
charged proton in its nucleus is felt by an incident electron when it comes 
into contact with a hydrogen atom. In comparison to the KN electronic cross 
section, which solely takes into account the scattering of electrons without 
taking into account the influence of a bonded nucleus, this attractive force 
increases the scatter probability, resulting in a greater KN cross section for 
hydrogen. 


Overall, the differences in the KN cross sections are caused by the interaction 
of atomic size, atomic number, and electron binding effects. Platinum has a 
bigger cross section than hydrogen and electrical interactions due to its bigger 
dimensions and higher atomic number. In addition, compared to electronic 
interactions, the effects of hydrogen's electron binding contribute to its 
greater cross section. 


It is found that the KN cross section in the relativistic case is greater than in 
the non-relativistic case when compared the KN cross sections in the 
relativistic and non-relativistic instances, as shown in figure 5. The scattering 
is governed by classical electromagnetic theory in the non-relativistic 
situation, where the electron energies are substantially lower than the speed 
of light (c), Assuming the amount of electron mass stays constant during the 
scattered process and ignoring relativistic effects, Thomson scattering is used 
to characterize the cross section for this case. On the other hand, relativistic 
effects become important in the case where the electron energy approach or 
surpass a large proportion of the speed of light. The KN formula, which 
integrates relativistic adjustments to the scatter amplitude and accounts for 
the fluctuation of electrons mass with energy, describes relativistic scattering. 


o--rm> Electronic KN Cross Section at a 
«w= Electronic KN Cross Section at a A 


oz Pt-Atomic KN Cross Section at a 
—Pt-Atomic KN Cross Section a, 
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Figure 5: KN cross section with energy ratio for relativistic and non- 
relativistic case 


Additionally, the electron's mass greatly increases in the relativistic zone as 
its kinetic energy approach the speed of light. The enhanced contact between 
an incoming electron and the intended atom caused by the higher mass has 
an impact on the scattering process. A higher scatter probability and, as a 
result, a greater KN cross section are produced by this increased contact. The 
KN formula considers the energy dependency of the scattering process while 
calculating the relativistic scattering amplitude. When an electron is struck at 
a greater energy, its wavelength approaches or even falls below the dimension 
of the scattering centers (atoms). Comparing to the non-relativistic scenario, 
this enables more effective scattering and increases the cross section. As a 
result, the relativistic situation has a greater KN cross section than the non- 
relativistic case due to the combination of the relativistic mass increases and 
the dependent on energy scattering amplitude. The increased cross section is 
due to relativistic phenomena that increase the scattering probability and 
enable stronger interactions among the incoming electron and the target atom. 


V. Conculsion 
In conclusion, the energy dependency of scattering processes can be used to 
explain the declining tendency of the KN cross section increasing electron 
energy. In comparing the cross sections of hydrogen and platinum, it is clear 
how important atomic size and number are; platinum has a higher cross 
section because of its bigger dimensions and greater atomic number. 
Additionally, the KN cross section remains constant in the high-energy area, 
indicating weaker and fewer energy-dependent interactions, possibly as a 
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result of relativistic processes. The KN cross section is found te be greater in 
the relativistic domain when compared the relativistic and non-relativistic 
situations. This is attributable to the relativistic mass rise and energy- 
dependent scattered amplitude. A larger range of atomic sizes, numbers, and 
substances should be investigated in order to better comprehend KN cross 
sections. Radiation therapy, medical imaging, and materials science can all 
benefit from the findings. Computer simulations are useful for understanding 
complex systems. 
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Abstract 


After three generations of quarks-leptons and Higgs boson have been 
confirmed, we research some possible development directions of particle 
physics. First, particles of three generations are extended. Second, the 
hadronic theory must precision quantification such as heavy flavor hadrons 
and the lifetime formula, etc. Third, various unifications of interactions are 
discussed. Fourth, some basic principles (Pauli Exclusion Principle, duality, 
uncertainty principle, etc.) are probably violated under some cases. Fifth, 
Santillli, Penrose, et al., proposed some new mathematical methods. Sixth, 
particle astrophysics, the negative matter as unified dark matter and dark 
energy, and the extensive quantum theory are searched. 


Key words: hadron, three generation, unification, basic principle, mathematics, 
astrophysics. 
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I. Introduction 
It is well-known that the three generations of quark-lepton in the standard 
model is: 
q=(2/3)e -(1/3)e -e 0 


u d e\ (vy, \ 

‘ é (1) 
wl | Y, 

t b T Vv 


So far these quarks-leptons and Higgs boson are confirmed. In this model the 
“ second and third generations of quark-lepton are the different excited states of 
the first generation. The standard model implies that third generation quarks 
are already final from bottom to top. 

After the basic elements and structure of the particle are determined, what 
is the further development direction of particle physics? In this paper, we 
research some possible directions. 

2. Extension of Three Generations 

The extension of three generations may be: 1. The most common method 
is to further confirm those new particles composed of the standard model, but 
these particles not yet being discovered. Different generations of quarks have 
a certain periodicity, most notably heavy quark charmonium and bottomonium 


7, =cc and Y= bb , which agree with Cornell potential [1]: 
mal a (2) 
r 


For c&, bb and i7, a,=0.35, 0.22 and 0.11. They become smaller with the 
bigger quark mass. Further, the mesons with c-quark are Dt =dc, D° =iic 
(l=1/2), and D* =c¥ (1=0). The mesons with b-quark are B* =ub , B® =db 
(l=1/2), and B =sb (1=0). Both are complete symmetry. The baryons with c- 
quark are (A‘,=*)=ude, **=uuc, L2=dde, Et =usc, B2=dsc, and Q2=ssc. 
Other E** =ucc, Bi =dec, Qt, =sec , QE, =ccc and so on. In 1990s, the 


heavy quark symmetry and heavy quark effective theory (HQET) are proposed 
[2-7], in which interactions between heavy quarks are mainly achieved by 
exchanging soft gluons and other light-flavor particles inside the heavy hadron. 
Usually phenomenological models are combined by some nonperturbative 
QCD. 

2. Find the gluon and whether it has three generations, or some excited 
states of the known first generation. 3. Whether W-Z have three generations, 
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or which are only excited states of W-Z. 4. Are there three generations of 
Higgs particles? 

It is known that K is weaker than z interactions. If it is generalized, D and 
B interactions should be weaker. Further, the three-generation interactions are 
getting weaker and weaker. If the interaction is the same, three quark-lepton 
exchange particles must be zero mass. 
3. Precision Quantification of Three Generations of Hadronic 
Theory 

The first quark-lepton and QED is very precise, which includes the 
anomalous magnetic moments of electron by the renormalization theory, and 
Lamb shift. Second and third quark-lepton should first be specifically 
renormalized. This corresponds probably to the Rosen-Ross mass formula of 
leptons-meson [8,9]: 


M=m,+—). (3) 
2a 
We proposed a modified accurate mass formula [8,10-12]: 
M = M, — AS + BLI(1 +1)-S?/2]. (4) 
For the J’? =1*/2 baryon octet, let M, =910.75MeV, A=222.04MeV and 


B=38.43MeV, so 
m(N) = 939.5725, m(A) =1115.6MeV, 


m(Z°) =1192.46, m(°) =1314.89MeV. 


For the J” =0™* meson octet, let A=0, 14,=549.4MeV and B= -207.22 MeV, 
so 


(5) 


m(1°) =134.96,m(K°) = 497.6, m(7) = 549.4MeV . (6) 
The neutral baryons and mesons agree completely within the range of error, 
and are unified M=m. 
According to the symmetry of s-c quarks in the same generation, the 
heavy flavor hadrons which made of u,d and c quarks should be also classified 


by SU(3) octet and decuplet. Such we derived the symmetrical mass formula 
[10-12]: 


M =M,+AC+ BUI +1)-C?*/2]. (7) 
Let M, =875, A=1453, B=84MeV, so m(N)=938, m(A‘)=2286 , 
m(z,,) = 2454 MeV, and we predicted m(&,,) =3676MeV and so on. In July 


ott 


2017 LHC observed a new and very charming baryon ©. =ucc, whose mass 
is 3621MeV [13]. Then m(S%,)=3621.4MeV, and J=1/2. It agrees with Eq.(7), 


whose error only is 1.4%. 
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For QCD, SU(3) has [7], 


m>. = Fm, +m,), (8) 


mM. =F(m, +m,), (9) 
m0 =F(m, +m). (10) 


By the experiments, (9)/(8) obtains m,/m, =11.507. (10)/(8) obtains 
m,/m,=11.70. Further, m,/m, =11.70/11.507=1.0168. But, (10)/(9) obtains 
m,/m, =0.01586. Both are complete different. Assume usually that 


m,/m, =2:1. It is that 4 parameters m 


nu? 


4v . 
m,,m, and ra determine 3 masses 


M _+,M,«,M,, and which is not agree with experimental data. 


Various masses of different quarks mass are very uncertain. And the 
quark mass change is irregular. Does it conform to the power M=m*, 
log M =klogm? 

Free quark not exist that can be considered in particle with the infinite 
deep potential well, which will only have several separate states E,, =n’E,. 


This should be related to the masses of the different quarks. However, only 
m(d):m (u)=4, and other does not conform. 

In the mass formula (4) the quantum number S corresponds to vibration, 
and | corresponds to rotation. C is also vibrations, both are basically 
symmetrical, A=228, and 1453. Quarks u-d only rotation l=1/2, s and c are I=0, 
do not form rotation. 

The Planck blackbody radiation formula and Bohr orbit are based on 
oscillators. The GMO mass formula corresponds to diatomic molecules with 
vibration-rotation. It corresponds to ud-s two bodies, and may develop to ud-c, 
etc. Assumed that rotation | is the same, while vibrations S and C are different. 
Further, S-C-B-T are all different vibrations, and probably correspond to 


n,n’,..n™ and so on. 


The heavy flavor hadrons include c, b, t quarks, but t quark is very 
unstable, and soon decays by the form t>b+W. But, its theory obtained only 
few quantitative results, and some introduce too many parameters. 

Based on the Y-Q and I-U symmetries between mass and lifetime on the 
general SU(3) theory, we obtained the lifetime formulas of hyperons and 
mesons [8,10-12]. Based on the new data, we proposed various lifetime 
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formulas of heavy flavor hadrons, which very agree with experiments. This is a 
new method on lifetime of hadrons described by quantum numbers, and can 
be unified for mass and lifetime [12]. But these still need to be precise and 
systematized. 

Present the mass formula is determined according to octet and decuplet. 


Whether there is a similar lifetime formulas? For octet p,n, Aj (I=0), =, (I=1) 
and =..(l=1/2) are a side of the SU(4). For decuplet three side masses are 
isometric or add one interaction term. m(Q™ )=1672.5, m(Q° )=2695.2, the 
mass difference is 1022.7. lsometry m(Q%, )=3717.9, m(Q** )=4740.6 (J=3/2). 
But, m( A** )=1232, m( x** )=2520, isometry m( 2%" )=3808, m( Q** )=5096 
(J=3/2). 

Moreover, m( A )=1232, m( =° )=2520, isometry m( Ei, )=3808, 


m(Q** )=5096 (J=3/2). Both m(Q:* ) are different, so we may determine that 


isometry is only approximate. 
Other aspects, various quantitative descriptions of decays and collides on 
particles still remains to be identified and developed. 
4. Unifications of Various Interactions 
In QCD, the function is [7] 
3 
& 2 5 
=— 11-=N,)+O(g’), 11 
B(g) ea 3 tg’) (11) 
where N, is the number of quark flavors. 1. If g>0, £ <0; if g<0, 6 >0; if g=0, 
£ =0. 2. For g>0, if N, <33/2, £ <0; if N, >33/2 which may not possible, £ >0; 
if N,=33/2, B=0.Usual f(g,) depends on the running coupling g,. 
127% 
a,(Q’) =—__———-____.. 
(33-2N,)In(Q’/A’) 
When Q?-»00 (higher energy and shorter range), a, =0, i.e., the 


asymptotically freedom. SU(3) is asymptotically free, which first is discovered 
by Gross and Wilczek [14] and Politzer [15] in 1973. 


B(g,) may be positive or negative or 0 [16-18]. A(g,)<0 (for O< g,<q), 
B(g,)=0 (when g,=q), and f(g,)>0 (for g,>q). From infrared attraction to 
ultra-violet repulsion [18], 


(12) 
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When energy Q? >~, the coupling constant of strong interaction a,(Q7)>0. 


Further, reduction of distance should be weak interaction with repulsive force. 
It can be inaccessibility and exclusion, are repulsion. 


We proposed that ((g,)>0 is strong interaction with big. g, and attraction, 
£(g,) <O is weak interaction with small g, and opposite repulsion, and 


&(g,) =0 with g, =q is the asymptotic freedom. The theory of strong 


interactions based on the gauge group [19-21]. This is also the unity of strong 
and weak interactions with short-range [19-21]. 
5. Possible Violations of Some Basic Principles 

So far the universality of Pauli exclusion principle (PEP) was queried 
some times. First in 1978-1980 Santilli proposed the rest of PEP, in particular, 
under strong interactions [22,23]. Based on some experiments and theories of 
particles at high energy, we investigated violation of PEP at high energy begin 
from 1984 [24-28], and researched some possible tests of violation of PEP 
(VIP) [25,27]. Recently, we researched another possibility on VIP at ultra-low 
energy, which is similar to superconductivity, superfluids and Bose-Einstein 
condensation (BEC), etc. Further, we proposed a possible mechanism of VIP: 
Cooper pairs extend to general fermion pairs, so they transform to bosons with 
VIP. From this we may predict some characters of this like-boson, as in nuclei 
and atoms, etc. Moreover, we researched some possible VIP in mathematics 
and physics. Based on the extensive quantum theory, its PEP may be violated 
[29]. At very low temperature two fermions can constitute a boson like the 
_ Cooper pair, and perform Bose-Einstein condensation. 

The ghost field has gone beyond the distinction of two kinds of quantum 
statistics, under what conditions does it hold? Is it at very high or at ultra-low 
energy? 


In 1995, the condensation numbers of ®’Rb and 7Li atoms may be high 
as 10° under this extreme condition [30]. In 1999 DeMarco and Jin cooled the 


potassium atomic gas with fermion characteristics to 10° K, so that the 
potassium atom pairs to realize quantum degeneracy to Fermi atomic gas [31]. 
Its quantum effects are different from Bose atomic gas, such as Fermi 
pressure, Pauli blocking and superfluid, etc. The interaction leads to the 
formation of Cooper pairs by Fermi atoms and the change of resonance - 
interaction to realize the phase transition from the supercurrent of BCS to BEC 
[32-34]. At ultra-low temperature fermion pairs can VIP. Contrarily, in 1960 
Girardeau proposed a gas model under the hard core boson limit, that is, 
Tonks-Girardeau gas, at this case the boson is confined to one-dimensional 
space and the repulsion is very strong. It is similar to PEP. In 2004 Paredes, et 
al., confirmed this case by cold atomic experiments [35]. This is a boson 
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similar to a fermion. So boson and fermion are symmetry and unification. In 
this case PEP has not played a role in the ultra-cold structure, and VIP may be 
tested. 

VIP Collaboration searched new experimental limit on VIP by electrons 
[36-38], and experimental tests of quantum mechanics on VIP and future 
perspectives [39]. Chakraborty, et al., discussed sufficient condition for the 
openness of a many-electron quantum system from the violation of a 
generalized PEP [40]. Abgrall, et al., researched new limits on bosonic dark 
matter, solar axions, PEP violation, and electron decay from the low-energy 
spectrum [41]. Shi, et al., searched experiments for VIP [42]. 

Based on the analysis of the logical structure [43], the duality, the wave- 
property is the basic principle of quantum mechanics. But, the wave property 
seems to have not been well tested under some conditions. We assume that 
wave should be corrected and developed, and for example, it is nonlinear 
wave, and a development of the linear superposition principle [8,44]. The 
quantitative restrict of wave property is mainly determined by 4=// p. For 


particles with high energy or large mass, the wavelength is particular small. 
The relativity is a range of high velocity and large momentum, while 
microphysics wave property exhibits easily only with low velocity and small 
momentum. We suggested that the basic characteristic of particles is the 
symmetry-statistics duality [8]. 

K. Shoulders (1991,1992,1999) found the charge cluster (Exotic Vacuum 
Objects, EVO): When 10 billions of electrons are compressed into an area of a 
millionth of a millimeter in diameter or less, they began to stick together. This 
not only violates the charge exclusion, but also fails PEP at high energy. It 
violates probably the charge and matter conservation laws. In 2000 Puthoff 
proposed a model attributed to the Casimir force produced by zero energy. 

Santilli proposed a unified form of possible generalizations of Heisenberg 
uncertainty principle for strong interactions [45]. Hadronic mechanics (HM) 
[46-48] has generalized the uncertainty relations in the only non-trivial way. 

In 2023 Santilli reviewed and updated that the insufficiencies of quantum 
mechanics in nuclear physics; the completion of quantum mechanics into the 
axiom-preserving, Lie-isotopic branch of hadronic mechanics for the invariant 
representation of extended protons and neutrons under potential and 
contact/non-potential interactions; the exact hadronic representation of all 
characteristics of the neutron in its synthesis from the proton and the electron 
at the non-relativistic and relativistic levels; the completions of Bell’s 
inequalities with ensuing iso-deterministic principle for strong interactions. 
Santilli then presented the apparent resolution of the historical objections 
against the reduction of all stable matter in the universe to protons and 
electrons and point out a number of open problems whose treatment is 
beyond the capabilities of quantum mechanics, such as: the cosmological 
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implications of the missing energy in the neutron synthesis, the prediction of 
negatively charged pseudo-protons, and the possible recycling of radioactive 
nuclear waste by nuclear power plants via their stimulated decay [49]. 

We researched uncertain relations of entropy, information and 
thermodynamics, and proposed the speed of light and some fundamental 
constants should be uncertain. Further, we searched its possible development 
in the entanglement, superstring, and the extensive quantum theory, various 
matters and different phases, etc. Moreover, the irreversibility and statistics 
are discussed, and we derived operators of entropy, etc [50]. New measuring 
technology and some developed theories already challenge the uncertainty 
principle, which may not be held under certain conditions, and must be 
modified and developed. 

6. Some New Mathematical Methods 

From 1978 Santilli discovered the new Iso-Mathematics, which includes 
isonumbers, isorepresentation [51-53], etc. In lso-Mathematics, isonumbers 
and genonumbers of dimensions 1, 2, 4, 8, and their isoduals and 
pseudoduals may extend to "hidden numbers” of dimension 3, 5, 6, 7 [52]. 
They are applied for many regions [53,54]. In 1998 Santilli shown that the 
objections against the Einstein-Podolsky-Rosen (EPR) argument are valid for 
point-like particles in vacuum (exterior dynamical systems), but the same 
objections are inapplicable (rather than being violated) for extended particles 
within hyperdense physical media (interior dynamical systems) because the 
latter systems appear to admit an identical classical counterpart when treated 
with the isotopic branch of hadronic mathematics and mechanics. Now Santilli 
reviewed, upgraded and specialized the basic mathematical, physical and 
chemical methods for the study of the EPR prediction that quantum mechanics 
is not a complete theory. This includes basic methods [55], apparent proof of 
the EPR argument [56], and examples and applications, in which the validity of 
the EPR final statement is the effect that the wavefunction of quantum 
mechanics does not provide a complete description of the physical reality. The 
axiom-preserving “completion” of the quantum mechanical wavefunction due 
to deep wave-overlapping when represented via isomathematics, and shown 
that it permits an otherwise impossible representation of the attractive force 
between identical electrons pairs in valence coupling, as well as the 
representation of al// characteristics of various physical and chemical systems 
existing in nature [57]. Moreover, Santill studied the classical determinism of 
EPR prediction by isomathematics [58]. 

It is known that field theory has scalar, vector, tensor, spinor, torsion and 
sO on. Penrose proposed twistor [59]. Penrose and MacCallum proposed 
twistor theory as an approach to the quantisation of fields and space-time [60]. 
Penrose searched nonlinear gravitons and curved twistor theory. Twistor 
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space is a parameter space-time with complex structure [61]. Penrose and 
Rindler discussed spinor and twistor methods in space-time geometry [62]. 

The space-dimension of mathematics and physics can be extended to 
higher n-dimensions geometry [63]. It may be Hilbert space and corresponding 
quantum mechanics, and superstring. In mathematics, physics and many 
regions the field theory is all a very important problem. When the space- 
dimension is extended to n, fractal and complex, and various number-systems, 
the field theory and its formulas may be correspondingly extended. In these 
cases, Gauss’s theorem and Stokes’s theorem, and corresponding extensions 
on gradient, divergence and curl are searched [64]. Further, they are 
combined each other, and form multiple combinations, such as the scalar- 
tensor fields, the scalar-spinor fields, the vector-spinor fields, etc. These fields 
can be applied to physics, biology, earthquakes and social science, etc. Field 
theory has been widely applied in many regions of natural and social sciences, 
and its any development will necessarily inspire and apply to more aspects 
[65]. 

In quantum mechanics the time-space operators of energy-momentum 
representation in quantum mechanics are [66]: 


xy ah Ly. (14) 

Op, 
Based on Eq.(48) we discussed the space-time operators and_ their 
generalization, and some applications of this method. In particular, the lifetime 
formulas of particles are obtained from the time equation, and they agree 
better with the experimental data. We proposed some operator equations of 
general relativity and special relativity. This is the simplest unifying quantum 
theory and general relativity, and corresponds to the extensive quantum theory, 
and may overcome the singularity problem in general relativity. It is the 
combination and unification on quantum mechanics and general relativity [67]. 

Further, based on the non-commutation 
[4,B] = AB- BA=n #0, (15) 

of matrix, group, tensor and so on, we propose the mathematical quantum 
theory. If 77 is imaginary number, it will correspond to the extensive quantum 
theory. If 7 is real number, it will be development of quantum theory. 
Moreover, 77 may be complex number, etc. We introduce a similar wave 
function and corresponding operators, various similar quantum results will be 
derived. Next, we discuss its physical meaning and various applications. Third, 
based the general matrix we research mathematics of unified gravitational and 
electromagnetic fields. Fourth, we discuss the space-time equations and the 
simplest unifying quantum theory and general relativity. Further, we can 
combine general discrete mathematics [68]. 
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The nonlinear approach of quantum mechanics is continuously an 
important direction. We propose mathematically the basic nonlinear operators 


: Ow 
Pu Se a) (16) 
wT : 
so Klein-Gordon equation and Dirac equations turn out to be 
(F’o+T-m)p=-J. (17) 
Dirac equations are: 
¥, (FO, +0 y + py =f. (18) 


We obtained the corresponding Heisenberg equation. Then the present 
applied superposition principle is developed to the general nonlinear form. The 
quantum commutation and anticommutation belong to F and I’. This theory 


may include the renormalization, which is the correction of Feynman rules of 
curved closed loops [44]. We think the interaction equations must be nonlinear. 
Many theories, models and phenomena are all nonlinear, for instance, soliton, 
nonabelian gauge field, and the bag model, etc. The superluminal entangled 
state, which relates the nonlocal quantum teleportation and nonlinearity, 
should be a new fifth interaction. Moreover, the nonlinear effects exist possibly 
for various interactions, for single particle, for high energy, and for small 
space-time, etc. The relations among nonlinear theory and electroweak unified 
theory, and QCD, and CP nonconservation, etc., are expounded. Finally, 
some known and possible tests are discussed [44]. 

We proposed that a trefoil knot of left hand and a trefoil knot of right hand 
is possibly similar to two spins s=1/2 and -1/2 of three quarks (nucleons p and 
n, etc) [69]. Two topological structures should differ in their energy levels in the 
magnetic field. Basic mesons of two quark form a line, and s=0. This model is 
first applied to proton and neutron, then it may be extended to other particles. 
7. Particle Astrophysics and Extensive Quantum Theory 

The particle astrophysics is a new modern intersecting science. Perkins 
published a book Particle Astrophysics [70]. In book some of the most 
perplexing cosmic and astrophysical phenomena are inextricably intertwined 
with the quantum world of elementary particles. Cosmologists have developed 
a “concordance model” that accounts for the observed global properties of the 
universe and the cosmic structures that were probably seeded by primordial 
quantum fluctuations. The standard model of particle physics accommodates 
all experimentally observed properties of elementary particles. But it 
completely fails to account for key elements of the concordance model—the 
dark matter that dominates the dynamics of galaxies and the dark energy that 
accelerates the expansion of the universe. Nor can it explain the cosmic 
preponderance of ordinary matter over antimatter, an asymmetry necessary 
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for our very existence. Well-motivated extensions of the standard model 
suggest novel particles as dark-matter candidates. 

The Dirac negative energy state corresponds to a negative matter with 
some new characteristics, which are mainly the gravitation each other, but the 
repulsion with all positive matter. Such the positive and negative matters are 
two regions of topological separation in general case, and the negative matter 
is invisible. It is the simplest candidate of dark matter, and is unified with dark 
energy, and can explain their some phenomena. We researched possible 
quantum theory of dark energy. The negative matter is related with inflation 
cosmos and many worlds, and discussed eight possible tests of this 
hypothesis, in particular, the season effect [71]. 

Based on Dirac negative energy, Einstein mass-energy relation and 
principle of equivalence, we propose the negative matter as the simplest 
model of unified dark matter and dark energy. All theories are known, only 
mass includes positive and negative. Because there is repulsion between 
positive matter and negative matter, so which is invisible dark matter, and 
repulsion as dark energy. It may explain many phenomena of dark matter and 
dark energy. We derived that the rotational velocity of galaxy is approximate 
constant. Assume that dark matter is completely the negative matter, so we 
may calculate an evolutional ratio between total matter and usual matter from 
1 to present 11.82 or 7.88. Further, the mechanism of inflation is origin of 
positive-negative matters created from nothing, whose expansion is 
exponential due to strong interactions at small microscopic scales. The 
negative matter as a candidate of dark matter and dark energy is not only the 
simplest, and is calculable and testable [72]. Generally, cosmology must be 
related to particles, in particular, for early Universe [73,74]. 

According to Feynman's idea [75] and based on a new form of the Titius- 
Bode law [76,77]: 


r, =an’, (19) 
we developed a similar theory with the Bohr atom model, and obtain the 
quantum constants =(aGM,)"” of the solar system and corresponding 
Schrédinger equation [76,77]: 
ra 
iH =-1PV'y +U-Ow. (20) 
Further, we proposed the extensive quantum theory in which the formulations 
are the same with the quantum mechanics and only quantum constant h and 
corresponding basic quantum elements are different [77,78]. 
8. Discussion 
In 2005, CERN experiments concluded that a new high-temperature and 
high-dense substance state has been produced. The quark-gluon plasma 
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(QGP) predicted by the original theory should be an approximate ideal gas 
composed of weakly coupled quark-gluons, but the result is that this matter 
exhibits extremely strong interaction, nearly perfect liquid (SQGP) [79], and 
may be described accordingly by the string theory of twin ecology [80]. This 
seems to violate the asymptotic freedom. ; 

Basis of quantum biology and quantum consciousness and quantum 
society [81-83] must be the extensive quantum theory [76-78], which includes 
the extensive quantum biology [84-86], and the extensive quantum sociology 
[87,88]. General drugs are always quantized. Only intravenous infusion is 
continuous, which corresponds to fluid mechanics. Its mathematical foundation 
and extension are mathematical quantum theory [68]. 
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Abstract 


In more than a dozen examples, the coordinate transformation known as the 
Special Relativity Theory is applied. The exact application of this trans- 
formation results in disagreement with generally known natural phenomena 
and with 5 experimentally verified laws of physics. At the same time, there 
is a disagreement with the first principle of Special Relativity Theory. 
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1 Introduction 


In 1905, Albert Einstein published a paper ”On the Electrodynamics of 
Moving Bodies” (” Zur Elektrodynamik bewegter Korper”) {1], which is now 
known as the Special Relativity Theory (SRT). Today, this theory is ac- 
cepted by the mainstream of physicists. Challenging this theory threatens 
the professional careers of physicists. 

The mainstream physicists believe that there is ample experimental ev- 
idence for the validity of SRT. Nevertheless they continue to invest huge 
amounts of time and effort to validate relativity (special and general), e.g. 
LIGO [2]. 

The ”twin paradox” is well known [3] and less well known is the ” Bell 
spaceship paradox” [4]. 

Special isorelativity (SIR) - [5] which has been developed for the dynam- 
ics of extended particles contains paradoxes as well (see in particular Sect. 
8 on). 


2 SRT 


First, let’s quote the necessary parts of the text from the translated original 
SRT [6]. 


p. 41: §2. On the Relativity of Lengths and Times 

”The following reflexions are based on the principle of relativity and on 
the principle of the constancy of velocity of light. These two principles we 
define as follows: 


1. The laws by which the states of physical systems undergo change are 
not affected, whether these changes of state be referred to one or the 
other of two systems of co-ordinates in uniform translatory motion. 


2. Any ray of light moves in the” stationary” system of co-ordinates with 
the determined velocity c, whether the ray is emitted by stationary or 
by a moving body.” 
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p. 48: ”so that the transformation equations which have been found 
become 


where 
BHldgiawie® 
p. 48: §4. Physical Meaning of the Equations Obtained in 
Respect to Moving Rigid Bodies and Moving Clocks 
” We envisage a rigid sphere of radius R, at rest relatively to the moving 
system k, and with its centre at the origin of co-ordinates of k. The equation 
of the surface of this sphere moving relatively to system K with velocity v 
is 
é? 4 nf? oh C = R? 


The equation of this surface expresed in 2,y,2 at time t=0 is 


(Jawa 


A rigid body which, measured in the state of rest, has the form of sphere, 
therefore has in a state of motion - viewed from stationary system - the 
form of an ellipsoid of revolution with axes 


R,/(1 — v?/c?), R,R. 


Thus, wheras the Y and Z dimensions of the sphere (and therefore of 
every rigid body of no matter what form) do not appear modified by the 
motion, the X dimension appears shortened in the ratio 1 : /1 — v?/c?, ice. 
the greater the value v, the greater the shortening.” 


y? +22 = R? 


p.49. "It is clear that the same results hold good of bodies at rest in the 
”stationary” system, viewed from a system in uniform motion.” 
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From the described transformations of the coordinate system, we prepare 
relations for the transformation of the distance | between two points and 
time interval T: 


le = Bly; ly = ly; le = 1,3 7 = Bt where B = 1/,/1 — v?/c? - (1) 


3 Alien Bob 


To describe the implications of using SRT transformations, we use the ex- 
ample of a second coordinate system moving at v=0.87c in the ecliptic plane 
of the Earth. Thus @ (or y) = 0.5. 

To illustrate, let us use the idea of an alien Bob moving towards (or 
away from) Earth at velocity v. It is irrelevant whether Bob has technology 
unknown to us or whether we use the idea of an alien Bob moving away 
from Earth at velocity v due to the expansion of the Universe. 


3.1 Earth rotation effect 


According to SRT, Bob sees, in accordance with [6] i.e. (1), a rotational 
ellipsoid with an axis ratio of 1:2. The orientation of the flattening of the 
ellipsoid depends on the direction in which Bob is moving. 

We know that the Earth rotates once every 24 hours, but the ellipsoid 
does not rotate, it is fixed in the direction of Bob’s movement. 


3.2 Elevator on the equator 


In areas on the equator, e.g. in Ecuador, everything moves up and down 
during the day. Like being in a huge elevator. It moves away from the 
center of the Earth and moves closer again. This is repeated twice a day. 

The size of the movement is remarkable. The height difference (the 
difference in the distance between the surface of the Earth and the center 
of the Earth) is more than 3,000km. 

The maximum speed of upward movement is 825km/h, i.e. the height 
of 93 Mount Everests per hour, that means the rise to the height of Mount 
Everest measured from sea level to the top in less than 39 seconds. 
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And then fall back down at the same speed. And this is repeated every 
12 hours. It is noteworthy that the time when the ascent and descent occur 
does not depend on anything on Earth, but depends on the direction of 
movement of the coordinate system, that is, the direction of movement of 
the alien Bob. 

It is also remarkable that the air pressure does not change during such 
movement, so there is no risk of mountain sickness like on Mount Everest. 
At the same time, the height difference is 360 times greater than the height 
of Mount Everest. 

It is remarkable that all rocks shrink and expand again even at great 
depths without mechanical damage. Reducing the height of a granite or 
limestone massif by half and then increasing it to its original height without 
cracking the granite or limestone is incredible. 


3.3. The orbit of the Earth around the Sun 


Alien Bob also tracks the Earth’s orbit around the Sun. It is not almost 
circular in his coordinate system as we know it, but elliptical. At first 
glance, this is nothing extraordinary, the planets move in elliptical orbits. 

It is noteworthy that the orientation of the ellipse, i.e. the direction 
of their axes, is not related to the movement of the Earth or something 
from the Solar System, but depends on the direction of movement of the 
coordinate system, i.e. the direction of movement of the alien Bob. 

Another oddity is the fact that the Sun is not in one of the foci of the 
ellipse, as it should be according to Kepler’s 1st law, but is in the center of 
the ellipse. 

The eccentricity of the ellipse e = a? — b? where a=lau, b=0.5au is 
e=0.87au. This is the distance of the focus of the ellipse from the center of 
the ellipse. Because lau & 150 mil. km. So the Sun is about 130 million 
km away from where it should be according to Kepler’s 1st law. 

Also, the area described by the Earth’s line segment is not constant. It 
differs by more than 20% 4 times a year, so the movement of the Earth does 
not coincide with Kepler’s 2nd law. 

Analysis of the measured speed and acceleration of the Earth’s elliptical 
orbit brings another surprise. 
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The speed of the Earth’s movement changes in a ratio of 1:2 2 times 
per year, i.e. 2x per one revolution of the Earth around the Sun. This is 2 
times more often than is measured for bodies observed from Earth moving 
in elliptical orbits. 

Also noteworthy is the centripetal acceleration calculated from the cur- 
vature of the Earth’s trajectory. The acceleration drops to 4of the maximum 
value twice a year, precisely at the point when the Earth is closest to the 
Sun. The maximum curvature of the Earth’s trajectory is when the Earth 
is farthest from the Sun. This contradicts at least Newton’s law of gravity. 

According to the law of gravity, the gravitational force should be 4 times 
greater when the Earth is closest to the Sun than when it is farthest from 
the Sun. The curvature of the trajectory should be proportional to the 
gravitational force. Thus, according to Newton’s laws, the curvature of the 
Earth’s trajectory should be greatest at the time when the Earth is closest 
to the Sun. 


3.4 The train effect 


For trains traveling on the east-west track in the equatorial region, two 
effects occur (2 times a day).The wagons and the people in the wagons are 
reduced to half their height. Twice a day, the distances between stations 
are halved and at the same time the tracks are shortened (including the 
shortening of the distances between track connectors and between sleepers). 

For trains traveling on the north-south line, the effect of changing the 
track gauge occurs twice a day. So the tracks will approach each other at 
half the distance. At the same time, the length of the sleepers and the 
length of the axles of the wagons and locomotives also change. So a normal 
train becomes a narrow-gauge train twice a day. 


3.5 Subway effect 


In the case of the subway (metro) line in the east-west direction, the height 
of the tunnel is reduced twice a day. At the same time, the height of the 
subway train, the height of the people in the subway, and the height of the 
seats in the carriages are also reduced. At the equator it is halved, at higher 
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latitudes less. In London or Paris it is only about a third. 

In the case of the metro line in the north-south direction, the gauge of 
the tracks changes twice a day, the width of the tunnel and also the width 
of the metro carriages also change at the same time. In them, the width of 
the seats located perpendicular to the direction of travel changes, and the 
width of the seats located along the walls of the wagons does not change. 

Passenger dimensions also change. Either nose size or shoulder width. 


3.6 Bridges effect 


Bridges are lowered and raised again twice a day. The height of the bridges 
is changing, including the height of the support columns. These change the 
more the higher they are. Although the support columns of different heights 
change differently, the bridge remains straight. 

If the bridge supports are made of concrete, they shorten and then 
lengthen without damaging the concrete. 

In addition, the bridges in the east-west direction change their length 
twice a day. Unlike bridges in the north-south direction, which change their 
width. In all cases, this is without damage to the material of the bridges 
and without external forces. 

E.g. Great concrete bridge in China [7], length 721m, length of the main 
field 445m, height above the valley 295m. Orientation east-west, position 
25°N. Twice a day, the concrete arch of the bridge lowers by approx. 135m 
and then rises again, while the ends of the arch come closer by more than 
220m. Although the bridge is made of concrete, it will not be damaged 
during such huge changes in dimensions. The deep valley over which the 
bridge leads does the same. 


3.7 Buildings effect 


Buildings also change their height twice a day, including lowering the ceilings 
and lowering the height of their residents. 

At the same time, buildings change their width or length twice a day. 
Windows in the south or north facade change their width and height, unlike 


- 38 - 


windows in the east or west facade/wall, which change their height and glass 
thickness. 

And that twice a day and without damaging the glass or frames. E.g. 
a window in a house in the equatorial area with dimensions of 2x2m in the 
south facade will decrease in height by 1m twice a day and narrow by 1m 
twice a day. And without damaging the glass or frames. 

E.g. the tallest building in the world, Burj Khalifa, 828m high, low- 
ers twice a day and rises again by more than 350m without any damage. 
The relatively smaller reduction is a result of the building’s off-equatorial 
location. 


3.8 Ship effect 


Ships change draft (and height) twice a day, including large seagoing ships. 
The height of the locks and the depth of the canals (e.g. Panama, Suez) are 
changing. During handling in the port, the behavior of ships changes twice 
a day. 

Twice a day, when the ship is turned by 90°, the ship shortens and widens 
or, on the contrary, lengthens and narrows. The draft does not change. 

2 times a day it will decrease (including reducing the draft) and it will 
not change either length or width when rotated. 

For example, the passenger ship Oasis of the Seas [8], which has a length 
of 360m, a width of 60m and a draft of 9m. 

This ship will shorten by 180m or narrow by 30m at the equator de- 
pending on the turn of the ship. The draft changes twice a day by 4.5m. 


3.9 Carousel effect 


People sit on the carousel and it spins. Twice a day, the situation occurs 
that the skull of each person on the carousel is deformed at each revolution, 
so that during one revolution of the carousel, it narrows twice and widens 
twice, and also shortens twice and lengthens twice. 

So, for example, on the carousel near the equator, during one rotation of 
the carousel, the skull, including the eyes, is halved twice, and the length of 
the skull, including the nose, is preserved. This alternates with a situation 
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where the width of the eyes and skull is normal, but the length of the skull 
is halved, the nose is halved, and the eye sockets are shallower and the 
eyeballs are flattened. 

This is replaced in 6 hours by a situation where the width and length 
of the passengers’ skull does not change when the carousel rotates, but is 
deformed from above, ie. it is half the height of the forehead, half the 
height of the eyes, half the teeth. 


3.10 The effect of high-voltage distribution systems 


Let’s observe the transmission routes of high-voltage power lines in the east- 
west direction. 

Twice a day the columns will decrease and increase again. Twice a day 
the distance between the pillars is shortened and increased again. 

When the columns get closer to each other, the wires will not sag more, 
but the opposite. The lowest point is higher than when the columns are fur- 
ther apart. The length and also the shape of the wires changes. When the 
columns are closest to each other, the cross-section of the conductors is cir- 
cular, when they are farthest from each other, the conductors are flattened 
with an elliptical cross-section. 

Let us recall that the resistance of a conductor is directly proportional to 
the length of the conductor and inversely proportional to the cross-section 
of the conductor, i.e. if the conductor is lengthened, its resistance will 
increase, if the cross-section is reduced, the resistance will also increase. 

In the case of high-voltage lines near the equator, the situation is as 
follows. When the distance between the columns is halved, the length of 
each conductor is halved and thus the resistance drops by 50%. 

In case of flattening of the conductor, i.e. in 6 hours after the shortening 
of the conductor, the cross-section will be reduced by half and thus the 
resistance will double. 

This means that the resistance of the conductor changes twice a day in 
a ratio of 1:4. 

The power loss in a conductor at constant current is directly proportional 
to the resistance of the conductor. So the power loss on the line conductors 
should change twice a day in a ratio of 1:4. 
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This ratio of resistances and therefore power loss depends on the speed 
of movement of the coordinate system, i.e. on the speed of movement of the 
alien Bob. The measuring instruments of the distribution network operator 
on Earth do not shew this. 


3.11 Glasses effect 


Two people in the equatorial region have glasses, beautiful glasses with 
round frames. Bob observes what is happening with their glasses. One of 
them faces west and the other faces north. 

For those who look to the west, the thickness of the glasses changes 
twice a day and thus the shape of the lens, which determines the optical 
parameters of the glasses (focal length / optical power - diopters). This 
alternates with a period when the shape of the frames, including the glass 
of the glasses, is deformed from circular to elliptical (recumbent ellipse). 
This is also reflected in the optical properties where the optical parameters 
of the glasses are different in different directions (cylindrical lens effects). 
Thus, the number of diopters of the glasses is reduced twice a day, and 
this alternates with the flattening of the frames and an increase in the 
optical power (increasing the number of diopters) depending on the direction 
(glasses with toric lenses compensating for astigmatism). 

The glasses behave differently for those who look to the north. The 
shape of the frames changes from a, lying ellipse to a standing ellipse. Four 
times a day, the frames of the glasses have a circular shape. 

The thickness of the lens glass does not change. The optical properties 
of the glasses also change. The axis of the cylindrical component rotates by 
90°(in steps four times a day). Cylindrical component in one direction in- 
creases, then decreases, axis changes, followed by growth and then decrease 
in the other axis and again direction change, increase and decrease... 

Like the lens of glasses, the lens of the eye and the eye deform. When 
the optical power (number of diopters) of the glasses decreases, the optical 
power of the lens of the eye and the eye itself also decreases. Optical effects 
add up. 

None of the users of the glasses subjectively noticed a change in the 
optical properties of the glasses. 
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3.12 Formula 1 races 


Cars driving on a race track have interesting behavior. This depends on 
direction as well as time (due to the rotation of the Earth). 

For example, let’s look at what Bob observes while watching the famous 
24 Hours of Le Mans race. 

During the race, the shape of the cars on the east-west sections of the 
track changes from short with tall elliptical wheels to long and low with flat 
elliptical wheels. 

On the north-south sections, low cars with flat wide wheels will change 
to narrow cars with narrow round wheels. The shape of the car and the 
wheels change when cornering because the direction of movement changes. 

Let’s take a closer look at the behavior of racing car wheels. That is, 
for spinning elliptical wheels. The car does not rattle as it should if the 
wheels were not circular. The ellipse does not rotate with the wheel. Only 
the wheel rotates and it deforms during rotation so that the flattening of 
the ellipse does not change direction. 

It’s easy to see if we look at the wheel valve. At maximum speeds, 
the wheels of the car rotate approximately 50 times per second. So the 
valve ” circulates” the center (bearing) of the wheel 50 times per second. 
The wheel disc and tire are deformed 100 times in 1 second so that the 
valve moves closer to the center of the wheel and moves away again at 
exactly twice the wheel rotation speed. When the car drives slower, the 
valve approaches less often, when faster, more often. At the same time, the 
disc of a racing wheel is made of very high-quality and strong materials, 
which require enormous forces for such large deformations. 

When Bob chooses a suitable spot on the Le Mans track, he can observe 
the effect of flattened (elliptical) wheels with a large flattening that occurs 
many times per second. And in the same place, the behavior changes during 
the race. 

When Bob chooses one car and follows it around the circuit, he notices 
that both the car and its wheels are changing. In addition, the transforma- 
tion also changes during the race. 

A short car with tall elliptical wheels turns into a long car with flat 
elliptical wheels and also a narrow car with narrow round wheels. 
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And since the Le Mans circuit lies at 48° latitude, these effects do not 
have the same magnitude (length ratios). 


3.13 Free fall effect 


Bob observes a, repeated free fall on the equator (e.g. in Ecuador) performed 
using a measuring device allowing the measurement of the free fall time from 
a height of 5m. From the measured values (fall length and fall time), Bob 
calculates the gravitational acceleration on Earth. They find that repeated 
free falls performed in the same location with the same equipment give 
different results. 

At the moment when the measuring device is farthest from the center of 
the Earth, it observes a free fall from a height of 5m, which lasts according 
to Bob’s clock Ar = 0.505s and has a gravitational acceleration of g = 
39.2m/s?. 

At the second moment, 6 hours after the first experiment, when the 
measuring device is closest to the center of the Earth, it finds that the fall 
from a height of 2.5m also takes Av = 0.505s. The resulting acceleration is 
g = 19.6m/s?. 

Thus, the gravitational acceleration changes in a ratio of 1:2 2 times a 
day and is the largest at the place that is farthest from the center of the 
Earth, and the smallest gravitational acceleration is at the place when the 
measuring device is closest to the center of the Earth. 


3.14 Pulley effect 


The use of the pulley has an interesting effect on the rope that passes over 
the pulley. On the equator, a situation occurs twice a day when the rope 
leading over the pulley lengthens or shortens during movement. 

When a rope is used that has marks on it that are equidistant from each 
other, and such a rope is guided over a pulley where it is almost horizontal on 
one side and vertical on the other (for example, when lowering and pulling 
a water bucket from a well), then the marks in the vertical part of the rope 
are closer together than in the horizontal part. When lowering, the marks 
on the rope that passed over the pulley get closer (the rope gets shorter). 
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On the contrary, when pulling out, the mark on the rope that passed 
over the pulley moves away. Both are without damage to the rope which 
is still taut. So the jack moves slower than the horizontal part of the rope 
moves, 


3.15 Tides 


Bob also observes the Moon and its orbit around the Earth. It is also 
elliptical with a semi-axis ratio of 1:2. The orientation of the flattening is 
the same as the flattening of the Earth and also the flattening of the Earth’s 
orbit around the Sun. 

According to Newton’s law of gravity, the tides should respect this ge- 
ometry of the motion of the Earth and the Moon. The maximum size of 
the tide occurs at the moment when the Moon is between the Sun and the 
Earth, that is, at the time of the new moon. Because the orbits of the Earth 
and the Moon are flattened, the distance of the Earth from the Sun and 
from the Moon at the time of new moon depends on the season. 

The amplitude of tides is inversely proportional to the cube of the dis- 
tance, according to Newton’s law of gravity. Therefore, the height of the 
highest tide should change in a ratio of 1:8 twice a year. A period of low 
tide should alternate with a period of up to 8 times larger tide after about 
3 months. 

Because the position of the Earth in its orbit around the Sun changes 
at the same time, the influence of the Earth’s flattening effect changes. At 
the time of the greatest approach, i.e. the greatest flattening, the height 
of the tide is also reduced on the equator to 1:2 due to the flattening. In 
Canada, the effect is smaller, so the height of the tide should be only 5.6 
times greater. 

Additionally, there is a length contraction effect that applies to a partial 
extent. So overall only 4x bigger. 

The highest measured tide is in Neah Bay, Nova Scotia, Canada [9], it 
has a size of about 11.5m. 

Due to the flattening of the Earth’s orbit and the Moon’s orbit, according 
to Newton’s law of gravity, a situation should occur twice a year that Bob 
will measure the height of the tide at about 11.5m and twice a year he 
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should measure the tide with a level difference of over 45m. 

It is noteworthy that the effect of buildings also occurs at the same time, 
i.e. the maximum tide according to Newton’s law of gravity is greater than 
a building 60m high. Because even the height of the building observed by 
Bob will decrease. 

This is not the case when Bob measures the sea level difference of Neah 
Bay so he measures the tide height in the range of 8.1 to 11.5m. At the 
same time, paradoxically, it will measure the maximum height of the tide at 
the time when the Moon and the Sun are farthest from the Earth and the 
minimum at the time when they are closest to the Earth. This is another 
contradiction to Newton’s law of gravity. 


3.16 The effect of the definition of the physical unit 
of distance 


In the history of physics, the distance unit meter has been defined in differ- 
ent physical ways: 

e 1791 - ten-millionth of the earth’s quadrant 

e 1799 - 1960 - the distance between the lines of the reference bar 

e 1960 -1983 --a multiple of the wavelength of krypton radiation 


e 1983 - fraction of the speed of light (distance in 1s). 


Let us analyze the properties of these definitions. 


3.16.1 Earth Quadrant 


The circumference of an ellipse with an eccentricity corresponding to Bob’s 
speed is 0.77 times smaller than the circumference of a circle. Thus, the 
length of a meter according to the definition from 1791 will be 0.77 less 
and therefore all distances measured by this shorter meter will be 1.3 times 
larger. 

As the Earth rotates, the length of the meridian changes. The meridian 
in the direction of Bob will be shortened, the perpendicular meridian will 
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remain unchanged. The length of the meridian will change twice a day, so 
the length of the distance unit will also change twice a day. And therefore 
the distances measured by this variable reference meter will also change 
twice a day. 


3.16.2 The distance of the lines of the standard 


The length of the standard will change depending on the direction and the 
time. Since the standard of the meter shortens as well as the dimensions 
of the objects measured by it, the numerical values of the distances come 
out unchanged. And therefore the length contraction by SRT will not be 
measurable. 


3.16.3 Wavelength of light 


There are two possible interpretations of the application of SRT transfor- 
mations to the wavelength of light. 

It is a distance and thus distance transformations are applied. This 
means that the wavelength depends on the direction of the beam. This will 
change with time due to the rotation of the Earth. So the light will change 
its color (wavelength) and the effects will be the same as in the case of lines 
on the standard. 

The second interpretation is that the wavelength of light is given by the 
relation \ = c/f where f is the frequency of the light, ie. f = 1/To where 
To is the period of oscillation. 

Time is transformed by SRT and ¢ is constant. Time transformation is 
independent of direction. 

In this case, Bob measures half the period of the oscillation and therefore 
twice the frequency. Therefore, the wavelength of light is halved regardless 
of direction. 

This means that the effect of length contraction in the direction of move- 
ment will numerically be zero, the numerical value of the distance will be 
the same, thus contradicting SRT. On the contrary, the distances measured 
in the perpendicular direction will be numerically doubled, i.e. also contrary 
to SRT. 
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3.16.4 The speed of light 


According to the SRT axiom, the speed of light is constant in all coordinate 
systems. 

Time transforms. Time transformation does not depend on the direc- 
tion of movement of the coordinate system. The distance transformation is 
dependent on the direction of movement. 

Light travels 300 thousand km in 1s in the direction perpendicular to 
Bob’s motion. Bob will measure a time of 0.5s, i.e. the length of the meter 
in Bob’s coordinate system will be 2m. This means that the distances 
measured with such a meter will be halved, i.e. contrary to SRT. 

Light travels 300 thousand km in 1s also in the direction of Bob’s move- 
ment. Bob measures a time of 0.5s, and at the same time measures a 
distance of 150 thousand km. So the length of the meter in Bob’s coordi- 
nate system will be 1m. So the distances measured with this meter will be 
in agreement with the SRT. 

The length of a unit of distance depends on the direction and therefore 
also on time due to the rotation of the Earth. 


3.16.5 Summary 


Thus, all definitions of meter used in history, except for one variant of 
the interpretation, have the physical effect of the time dependence of the 
reference. That is, the time dependence of all measured distances using a 
physical unit defined in this way. 

The only exception is the interpretation based on the frequency of the 
reference light source. In this case, the length of the reference is independent 
of the direction and therefore of the rotation of the Earth. In this case, 
the length of the reference meter depends on the speed of movement of 
the coordinate system. In this case, the numerical values of the measured 
distances will contradict the SRT regardless of the direction. 


-47 - 


4 Alien Alice 


The alien Alice, unlike Bob, moves at a speed of v = 0.44 c, i.e. roughly 
half the speed of Bob. It also moves in the plane of the Earth’s ecliptic. For 
Alice, 8 or y = 0.9. 

In Alice’s coordinate system, all the described effects are visible as in 
Bob’s coordinate system, except that the numerical values are significantly 
different. 

For example, the maximum speed of upward movement due to the 
Earth’s rotation is 165km/h at the equator. That’s 5x slower than Bob 
observes. 

The eccentricity of the Earth’s orbit is 0.14 and thus 6 times smaller. 
The ratio of resistances of distribution network conductors is only 4:5, i.e. 
roughly 13 times smaller. A concrete bridge in China lowers ”only” by 
about 30m, the passenger ship Oasis of the Seas changes its draft by less 
than a meter. In both cases, 5 times less. 

Gravitational acceleration calculated from free fall experiments from 
data measured by Alice is 3.2x smaller than from data measured by Bob. 

When Alice communicates with Bob and they share the results of their 
observations, they discover that they are observing the same effects but the 
measured values are many times different, i.e. as if different physical laws 
work for different observers on Earth. 

Therefore, together they compile a table of contradictions of observed 
effects with experimentally verified physical laws. They find that the fol- 
lowing laws of physics are violated: 


1. 1st and 2nd Kerpler’s law 

2. Newton’s law of gravity 

3. Laws of strength and elasticity, including Hooke’s law 
4, Ohm’s law 

5. Laws of classical (geometric) optics 


This means that the 1st SRT principle is not fulfilled 
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5 Conclusion 


SRT coordinate transformations violate the first principle of SRT. SRT con- 
tradicts at least five experimentally verified laws of physics and practical 
experience of earthlings. 
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Abstract 

Once a number of fundamental thermodynamic notions associated to a reversible, 
ideal-gas phase chemical reaction were determined, among them the equation for its 
path, the equation for its driving force, and via this last one, the equation quantifying 
its work output at any given point along its unfolding, we proceeded to study in detail 
the occurrence of this reaction in a van’t Hoff equilibrium box. The results obtained 
not only evinced the similarity existing between the way work is produced in both, 
the van’t Hoff equilibrium box and Carnot’s heat engine, but also the fact that like 
the latter, the former is also a cyclical process. The fundamental difference between 
the two being that while in Carnot’s engine work is produced out of the heat 
transferred from a hotter to a colder heat reservoir, in van’t Hoff’s work comes out 
of the transfer of chemical matter between two different chemical states described 
by chemical potentials. 
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1. Introduction: The system of interest and its associated thermodynamic 
data, 

The arguments to be here presented will be centered on the following 
thermodynamically-reversible, stoichiometric, isothermal, and _ isobaric 
chemical reaction taking place in the ideal-gas phase: 


1 Butane (ideal gas, 1 bar, 298.15 K) 
= 1 Isobutane (ideal gas,1 bar,298.15K) (1) 


For the purposes of this discussion reversible is any process which in 
concatenation with its inverse restores the condition existing before any change 
had taken place, the so-called initial condition. 


In the thermodynamic analysis of reaction (1) performed in the first six sections 
of this paper, the reaction system, namely pure butane at the initial state, pure 
isobutane at complete conversion, and any intermediate reaction mixture 
constitutes the system. The surroundings are constituted by a heat bath, in charge 
of maintaining the system’s isothermal condition, and a mechanical or work 
reservoir acting as depository or supplier of the work either produced or 
demanded by the reaction system. These three bodies constitute the universe of 
the process, Serves to note that on reason of it taking place without change in 
total mole number, isothermal and isobaric reaction (1) is also isochoric, that is, 
it also evolves at constant volume. The relevant thermodynamic properties for 


this reaction are given in Tables 1, and 2. 
= 
“17.1141 | 


“eine | 03060 
-134.516 -0.3812 | 


Table 1 Relevant standard thermodynamic data for butane and isobutane at 
298.15 K (Stull et al., 1969, pp. 209, 245). AS? stands for the standard molar 
entropy of formation in kJ K~* mol, with AH? and AG, respectively 
representing the standard molar enthalpies, and standard molar Gibbs energies 
of formation, both in kJ mol~?. 


-20.8612 
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Reaction AH° AS° AG° 


kJ mol~* kJ K~2 mol~1 kJ mol~+ 


butane — isobutane | -8.3684 -0.0155 -3.7471 


Table 2. Standard enthalpy, entropy, and Gibbs energy changes for reaction (1) 
per mole of reaction. These values were obtained as the difference of the 
corresponding enthalpies, entropies, and Gibbs energies of formation of 
isobutane minus those of butane given in Table 1. 


2. The chemical potentials and the chemical force 
The Gibbs energy change for a homogeneous, closed system capable of 
undergoing changes in composition, can be written as follows: 

dG =VdP —SdaT —dW (2) 
The previous equation can be obtained from the combination of the defining 
equations for the internal energy dU = TdS — PdV — dW; enthalpy dH = 
dU + PaV + VdP; and Gibbs energy dG = dH — TdS — SdT. 


In eq. (2) dW represents work different from the work of displacement PdV, 
that is, different from: 


... the work done by the system in displacing its environment at the 
steady pressure P... (dW) is the chemical work ... a form of work which 
can be done by a system in the absence of any change of volume and due 
to its change of composition ...and which can be attained by use of the 
van’t Hoff equilibrium box or the galvanic cell.” (Denbigh, pp. 68, 81) 

If the evolution of the system is isobaric and isothermal, then eq. (2) takes the 
following form: 

-dG=dw (3) 
It is on reason of the state function nature of the Gibbs energy that the previous 
equation applies to both, the reversible as well as the irreversible transit of the 
system between identical initial and final states: 


—dG = dWrey = AWrost (4) 
Or, equivalently 
-—AG = Wrey = Wiost (5) 
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Equations (4) and (5) make clear the fact that for isothermal and isobaric 
processes the Gibbs energy change represents the work transferred to the 
mechanical reservoir along the reversible path (W,e,), aS well as the dissipated 
work-producing potential (Wj0;;) transferred in the form of heat to the heat 
reservoir along the irreversible path: 


The change in G is, of course, the same for the same initial and final 
states of the reaction system, whether the process is conducted reversibly 
or not. It is the output of work which varies ... (Denbigh, p. 76) 


A system undergoing change can perform maximum work only when the 
change is carried out reversibly. If the process is not completely 
reversible, the amount of work obtainable is always below the maximum, 
the difference appearing as heat. (Maron & Prutton, p. 189) 


If the system under consideration turns out to be a homogeneous mixture of two 
components, a and #, then the following equation, making explicit the 
dependence of its Gibbs energy from temperature, pressure, and composition, 
also applies: 


G=G(T,P,neng) (6) 


In what follows a is to be an alternative descriptor of butane, and £ of isobutane, 
as they appear in reaction (1). 


The total differential of eq. (6) has been written below: 
dG = (0G/AT)pn,dT + (0G/OP) rn, 4P + >.@6 Jandrprjdn (7) 
i 


The usual provisos apply regarding the n; and n,; constraints appearing in the 
partial derivatives of the previous equation. 


For a constant temperature and pressure evolution of the reaction system, the 
previous equation reduces to: 


dG =) (0G/An)rpnjam (8) 
t 


The comparison of equations (8) and (3) allows us to write: 


dG =-dW =) @G/An)renjam (9) 
t 
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The partial derivatives appearing on the right-hand side of equations (8) and (9), 
of capital importance in the thermodynamic theory of spontaneity, chemical, and 
phase equilibrium, are also represented in the following simpler manner: 


Hi =(0G/ONj)prn, — (10) 


This intensive quantity ,;, with units of energy/mole, is the chemical potential 
of species i in the reaction mixture. As such, its value is dependent on the 
pressure, temperature, and composition of the said mixture. 


The chemical potential of component i can be interpreted as the amount by 
which the capacity of the reaction mixture for doing work (other than work of 
expansion) increases per mole of i added at constant temperature, pressure, and 
composition. The following example provided by Denbigh might serve to make 
its meaning clearer: 


...the chemical potential of copper sulphate in its acqueous solution in a 
Daniel cell is equal to the increased capacity of the cell to provide 
electrical energy, per mole of copper sulphate added (the addition 
actually being an infinitesimal). (Denbigh, p. 79) 
Even if the ‘infinitesimal addition’ referred above by Denbigh is a way of 
satisfying the constant composition constraint implicit in eq. (10), a more 
sensible way to this end would have consisted in adding a small fraction of a 
mole of copper sulphate to a finite quantity of solution followed by the 
calculation of the proportionate increase per mol. “If (this) is done, for strict 
accuracy one must pass to the limit of (AG;/4n;), as dn; approaches zero for 
the substance in question. The limit of this ratio, however, is exactly uj: wy = 
jim G;/4n,)” (Aston &Fritz, p. 89; The quote includes minor changes in 


notation.) 


In terms of the chemical potential, the following equivalent expressions for 
equations (8) and (9) will be written: 


dG = y mdm (aa) 
i 
dw = -) udm (12) 
i 


Let us now bring here the thermodynamic variable known as the degree of 
advancement, or reaction coordinate ¢, defined as follows (Moore, pp. 283-284): 
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An 
=— mol (13) 
Vi 
In the previous equation An; represents the change in number of moles 
experienced by chemical species i at the point at which ¢ is being evaluated, 
with v; representing its unitless stoichiometric coefficient, negative for 
reactants, and positive for products. For reaction (1) it adopts the following 
forms: 

An; Ng-1 ng—O 

f= = (14) 

Vi —1 1 

From the previous expressions we get: 
Ng =1-E, ng=s (15) 
The fact that the total number of moles is equal to _ one: 
Ng tng =1-F+F=1 (16) 

allows us to realize that the mole fractions (¥;) of these chemical species in 
the reaction mixture, as well as their respective partial pressures (p;) in it, 
this on reason of the total pressure being equal to 1 bar, are numerically 
identical to their respective number of moles, as given by eq. (15), i.e., 


Nq =Va=Pa=1-F, nmge=Yp=pp=s (17) 
The corresponding differential expression for eq. (13) is given below: 
dn; = v,dé (18) 
The substitution of eq. (18) in eq. (12) leads to: 


aw = (->. vat) dé (19) 


L 
In terms of the definition of work, thermodynamic or otherwise, as the product 
of force times displacement (Klotz & Rosenberg, p. 45), we can express dW as 
follows: 
dW = chemical force x reaction’s displacement = f d& (20) 


A simple comparison between equations (19) and (20) allows us to identify the 
chemical force as the expression between parenthesis in eq. (19), i.e., 


f=-Yvm 21) 


i 
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When the properly signed stoichiometric coefficients vg = —1 and vg = 1 for 
reaction (1) are substituted in eq. (21), we get: 


dW =(1u_—1mg)d& = (22) 
The previous equation permits writing the following expression for the chemical 
force driving reaction 1 Butane > 1 Isobutane: 
f = 1g — lug (23) 


The following equation, coming out of the combination of equations (20) and 
(22), can now be written for reaction (1): 


dW = (1g —1ug) dé = fag (24) 
The fact that dé is, by definition, a measure of the change in the number of moles 
of reactants and products allows us to bring here the interpretation provided by 
Brenstead of the dW term appearing in the previous equations, as “... the work 


associated to the transport of chemical matter between two different chemical 
states described by chemical potentials.” (Gill, p, 507) 


The following statement of Maron & Prutton (p. 189, italics are ours), asserting 
that: 
Work results only when the tendency of systems to attain equilibrium is 
harnessed in some way. From a system in equilibrium no work can be 
obtained, but a system in its way to equilibrium may be made to yield 
useful work. 


will be given here the following representation: 


dw 
Ce) 2g 605) 


eq 
Combination of equation (25) with eq. (24) leads us to: 
dw 
(=) = lie eq = 1B ,eq = teq = 0 (26) 
dé ee 


The previous result allows us to trace this (€W/dé) eq = 0 characteristic of a 
reaction in equilibrium to an even more fundamental condition: the fact that at 
this point the chemical force f—and with it the difference in chemical potentials 
it represents— is equal to zero, i.e., 


Teq = 1taeq - Lup eq =0 (27) 
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Once the chemical force has been re-expressed as a function of &, then the work 
associated, say, to the transit (0 — €) of the reaction system, can be evaluated 
via the following integral: 


é 
Work = | fede (28) 
0 


3. Thermodynamic work 


Work, says J. C. Maxwell (p. 54), “...is the act of producing a change of 
configuration in a system in opposition to a force which resists that change.” 


Regarding the work of expansion of a gas, the fact that “... for a given expansion 
the work performed depends upon the external pressure resisting the 
expansion...,” makes it necessary to distinguish the (internal) pressure of the 
gas, P, from the external one, P,,.. “In order that an expansion may occur, the 
internal pressure, must, of course, be greater than the external, but this difference 
may be made as small as is desired, and as a limit we may consider the case in 
which the internal pressure is equal to the external pressure or differs from it by 
a negligible amount.” The work produced under this condition is the reversible 
work. (Quotes from Pitzer and Brewer, p.35) 


In general, while the work of expansion can be represented as W = f P,,.dV; the 
reversible work, on its part, takes the following form W,e, = f P(V)dV. 


The previous statements require the following clarification: 


While thermodynamics cannot do away with the infinitesimal difference between 
driving and resisting forces in reversible processes —making them equal makes 
the process impossible; mathematics certainly can, this on reason of the fact that 
the numerical results obtained with one or the other are for all practical 
purposes, identical. 


This seems the proper place state that reversible processes are ideal, that is, 
practically impossible processes, and this is on reason of the impossibility of 
achieving an infinitesimal difference between driving and resisting forces. It is 
not uncommon to find the impossibility of reversible processes also defined in 
the following manner: “Strictly speaking, all reversible processes are impossible 
in nature, as they would require an infinite time for their accomplishment.” 
(Maron & Prutton, p. 110) Even if practically unrealizable, the fact that they 
represent the efficient limit to natural processes, makes of them an indispensable 
tool for the theoretical analysis of thermodynamic processes. 


Now, the following words from Denbigh (p.78): 


257 = 


The chemical potential has an important function analogous to 
temperature and pressure. A temperature difference determines the 
tendency of heat to pass from one body to another and a pressure 
difference determines the tendency towards bodily movement... a 
difference of chemical potential may be regarded as the cause of a 
chemical reaction or of the tendency of a substance to diffuse from one 
phase to another. The chemical potential is thus a kind of ‘chemical 
pressure’ and is an intensive property of a system, like the temperature 
and pressure themselves, 


allows us to make a simile (with all the appropriate reservations) between the 
chemical force given above by eq. (21), and the internal pressure of the gas. The 
internal qualifier comes on reason of the fact that the chemical force finds 
definition in terms of the chemical potential of the substances constituting the 
reaction system; it is thus, as Denbigh notes, a property of the system, just like 
P isa property of the gas. 


Just like the equation W,ey = f P(V)dV is used under the assumption that a 
pressure infinitesimally smaller than that of the gas is resisting the expansion, 
the same can be done with equation Wye, = { f (&)dé to evaluate the chemical 
work. This procedure de facto exempts us to provide any details as to the nature 
of this resisting force, or about the mode in which it performs its opposing or 
resisting role. Further comments on this issue in the discussion section. All we 
need to know to perform these integrations is knowledge of the functional 
relation connecting P and V in the former, and f with é in the latter. The 
interesting thing to be noted here is that even if f is the explanation for the 
occurrence of the reaction, and fdé is the quantifier of its reversible work output, 
this pair plays, however, no role in the actual production of this work. As will 
be seen below, for chemical reactions taking place in the van’ Hoff equilibrium 
box, this task, the actual production of work is carried on by the pair PdV. A 
similar thing happens in heat engines. While it is a temperature difference the 
one driving the flow of heat from the hot to the cold reservoir, work appears in 
the surroundings via the net result of two isothermal and reversible volume 
changes: one expansion, and one compression. 


Even, if the work defined by the integral We, = f f(&)dé is, as made explicit 
by the ‘rev’ subscript, reversible work, it also quantifies, for an irreversible 
transit between the same initial and final states, the work lost in it. Being this so, 
we will abstain from attaching any superscript to W. 
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4. The functional connection between chemical force and reaction extent or 
degree of advancement in reaction (1) 

According to Denbigh (p. 115): 


... a gaseous mixture (is) said to be perfect if the chemical potential of 
each of its components is given by the following relation, in which pu? is 
a function of temperature only, 
Mi =u? + RT InP + RT Iny; 
Where P is the total gas pressure and y; is the mole fraction of component 
i. Now, since ? is independent of composition it retains the same value 
when y; is brought up to unity. It is thus precisely the ... value of the 
Gibbs free energy per mole of the gas i in its pure state at unit pressure. 
This equation ... can be put in a more compact form by means of the 
partial pressure p;. Thus 
Mi = up + RT np; (29) 


The combination of the previous equation with the relations established by eq. 
(17), allows us to write the chemical potentials of butane and isobutane, and the 
chemical force by them defined, as they all pertain to reaction (1), in the 
following manner: 


Me =e +RTIn(1-§) (30) 
Mg =egtRTInG) (1) 
f =1[u$+RTInG -€] -1[ugt+RTIn()] (32) 


After rearrangement and simplification, the previous equation becomes: 
: 
33 
faz) (8) 


Further replacement of Ay° for the standard Gibbs energy change of reaction 
(1), produces: 


f =—Ap? - RT In( 


—f =AG° + eTin(=+5) (34) 


In which case, and in accord with equations (4), (24), and (34), we can write the 
following equation: 


dG = —-dW = —fdé = [ace + RT In (5) dz (35) 
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Which can also be written as follows: 


dG dw : ¢ 
= a4 +RTIn(7+4) (36) 


Applying the equilibrium condition stated above in eq. (25), which given eq. (4) 
can also be expressed as Ca eq = 0, we obtain the well-known equation: 


AG? = —RT In on =—-RTInK (37) 


1-—¢ 
While @ eq = 0 expresses, on one side, that at the state of equilibrium, as the 


previous quote of Maron and Prutton states, no more work can be obtained from 
a reversible chemical reaction, nor dissipated by an irreversible one; the 


Co eq = 0 condition identifies, on the other side, the state of equilibrium as the 
minimum of the graph G — vs. —é: 
The criterion of equilibrium for a system of prescribed temperature and 
pressure is that G has reached its minimum possible value. (Denbigh, p. 
83) 
5. Numerical values for the chemical potentials and the chemical force for 
reaction (1) 


With R = 8.314 JK~1mol? and the respective standard Gibbs energies of 
formation of butane, and isobutane at 298.15 K, namely 
—17,114.1 Jmol™1 and — 20,861.2 J mol}, we will have equations (30), (31), 
and (32) taking the following forms: 


Mg = —17,114.1 + 2478.8In(1 —&) Jmol™+ — (38) 
Mg = —20,861.2 + 2478.81In(G) Jmol™* = (39) 
f = 1ptg — 1g = 3747.1 + 2478.8 1n (Fs) Jmol? = (40) 
The evolution of these magnitudes along the course of reaction (1) can be read 
from Table 3, constructed by the simple expedient of substituting a few selected 


values of & into each of the three equations given immediately above, and 
performing the indicated operations. 


In column 1 of Table 3, €¢ = 0.82 was calculated by solving eq. (37) once 
AG° = —3747.1 J mol~1, and RT = 2478.8 J mol~+ were substituted in it. 
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The combination of equations (28) and (34) produces the following equation for 
the work associated to reaction (1): 


W = | —[AG° + RT nga) dé (41) 

Performance of the previous indefinite integral produces: 

W =-€ AG? —RT[En(Z) + (1-H nA -H]—RTHC 
The constant of integration C in the previous equation, evaluated with the known 
fact that at complete conversion, € = 1, it is true that W = —AG°, turns out to 
be equal to C = RT, in which case, we can write it in the following manner: 

W =-£AG°-RT[EM(E)+A-HmMaA-H) (42) 
The substitution in eq. (42) of the thermodynamic magnitudes there involved 
produce, in turn, the following equation: 

W = 3747.1 € — 2478.8 [EIn(é) + (1 — €) nd - $)] (43) 

It was via the simple expedient of substituting in eq. (43) selected values for ¢, 
that the numbers for W appearing in columns 5 and 6 of Table 3 were obtained. 
6. About the results consigned in Table 3 


The results produced by equations (38), (39), and (40), for the chemical 
potentials of butane, isobutane, and for their difference, that is for the chemical 
force, at various degrees of advancement of reaction (1), are those respectively 
contained by columns 2, 3, and 4. In accord with their respective dependence of 
the concentration of butane and isobutane, we see fg decreasing, and pg 
increasing as the reaction progresses, up to that point, the state of chemical 
equilibrium, at which no more difference exists between them. In accord with 
these considerations, we see the chemical force transiting from an arbitrarily 
large magnitude at € = 0, to a value of zero at equilibrium. The work calculated 
with the previous equation, identical at every & with the negative of the 
corresponding Gibbs energy change, quantifies, as explained above, the work 
gained, that is the work available in the work reservoir at any given ¢ along the ° 
reversible path, as well as the work lost at same degree of advancement, along 
the irreversible path. The numbers in column 5 corresponding to post- 
equilibrium é’s, applying only to the reversible path, are smaller than those 
corresponding to gg = 0.82. The decrease obeys to the fact that work must be 
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removed from the work reservoir in order to be spent in forcing the system to go 
beyond its equilibrium condition, that is, into the region of conversions of 
isobutane larger than &,, = 0.82. Concomitant to the decrease in the work 
available in the work reservoir, we find, on reason of this work being done in 
the system, an increase in its Gibbs energy. For example, the 63 J decrease in 
the content of the work reservoir, from 4241 J to 4178 J, spent in pushing the 
system along the transition (¢,g = 0.82) — (¢ = 0.9), appear as an increase in 
the system’s Gibbs energy, from —4241J to —4178J J along the same 
transition. 


§ —1ug | —1ug Pad W =-AG W; —W, | 
ae OF 


17114.1 


0.10 17375.3 | 26568.8 | 9193.5 1180 1180 


oe ama fama esa [as 
Ce 


Table 3. Column 1 in mol, Columns 2 through 4 in Joules/mol, and columns 5 
and 6, in joules. The graphs of these numbers are shown in Figures 2 and 3. 
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That the numbers in column 5 for § > &eg do not apply to the irreversible path 
obeys to the fact that at no point along this path, including the state of chemical 
equilibrium, we find any work available in the work reservoir. 


SSS 
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Chemical potentials , kJ/mol 


Figure 2 The chemical potentials of butane (a) and isobutane (A) vs. é. 


Chemical force , f kJ/mol 


Figure 3 The chemical force f is positive for §.g < 0, zero at feg = 0.82, and 
negative for eg > 0 
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It is this fact, the absence of available work in this irreversible universe, what 
makes any transit beyond the equilibrium condition, impossible. To make it 
happen we must re-define our universe to accommodate in it a new work- 
supplying body to carry on this task. 


That the chemical force f is, as shown by equation (36), the instantaneous rate 
at which the work-producing potential of the reaction, as measured by AG, is 
being transformed into actual work along the reversible path, or dissipated into 
heat along the irreversible one, explains the heading of column 4. 


The diminishing amounts of work produced (dissipated) as the reaction proceeds 
are the matter of the numbers in column 6. Thus, while the contribution of the 
transit from (§ = 0) > (€ = 0.1) amounts to 1180 J, the one from (€ = 0.1) > 
(€ = 0.2) reduces to 810 J; the one from (€ = 0.8) > (eq) only amounts to 3 
J. This trend is a consequence of the fact that the value of the chemical force 
decreases as the reaction proceeds. 


Let us now address the issue of how is it that the work producing potential of 
reaction (1) becomes gained work, that is, work available to us in a work 
reservoir. 


7. A thermodynamically-reversible chemical reaction in the ideal gas-phase 


7.1 The isomerization of butane into isobutane in the van’t Hoff 
Equilibrium Box 


The present section is a modified version of the one appearing in a previous 
paper by this author on a related subject (Iftiguez, 2014) 


The detailed description of the mechanism through which the efficient 
conversion of the work producing potential of a chemical reaction into actual 
work takes place, will be done in what follows by studying the occurrence af 
equilibrium of reaction (1). The device in which this reaction is to take place, 
known as the van’ Hoff equilibrium box (£B), has been represented in Figure 1. 


The operation of the EB for different types of chemical reactions, namely 
homogeneous reactions in the gas phase, and heterogeneous solid-gas reactions, 
has been described with various degrees of detail by Kauzmann, pp. 291-297; 
Denbigh, pp. 72-74; Smith and Van Ness, pp. 409-410; Fermi, pp: 101-106, and 
Pauli, pp. 54-56. among others. 


ee 


The salient features of a reaction occurring at equilibrium in a van’t Hoff 
equilibrium box, are the following: 


1) It will take place without the production of work; conclusion which in terms 
of our previous discussions, finds explanation in the fact that at equilibrium it is 
true that feg = 0, and if so, also true that dWeg = fegd§ = 0. For such a 
reaction, the expression of the first law, AU = Q — W, reduces to AU = Q, 
which means that all the energy released by the reaction ends up as heat in its 
heat bath. The bath will later exchange some of this energy with a number of 
isothermal and reversible volume changes involving reactants and products, for 
it to be transformed, via the action of the pair PV, into the thermodynamically 
predicted maximum possible amount of work. 


2) The final result of the operation to be described of the equilibrium box, takes 
the form of a standard reaction as in it the reactants remaining and the products 
produced appear in their pure states, that is separate from one another, each at 
the standard pressure of 1 bar at the temperature of the reaction. 


The thermodynamic analysis starts with the reaction box, designated as Cp in 
Figure 1, containing a vast amount of reaction (1) in chemical equilibrium at 
298.15 K and 1 bar total pressure. It is to this equilibrium mixture that 7 moles, 
0 <n 1, out of the initial mole of butane available for reaction (1), will be 
added, and n moles of isobutane removed. It is in this vast amount of reaction 
mixture in equilibrium where we will see the said n moles of butane (a) 
introduced in it, isomerized into n moles of isobutane (6). The (1 — n) moles of 
a remaining and the n moles of £ produced will be found at the end, separate of 
one another, and each of them at lbar and temperature T. 


The result of the four processes which described below, constitute van’t Hoff’s 
box mode of operation, can be represented as follows: 
1a (g,1 bar,T) > [((1 —n)a@ (g,1 bar, T)] + [nB (g,1bar,T)] (44) 

The brackets appearing in the right-hand side of the previous equation 
emphasize that these chemical species are found at the end of said four 
processes, separate from one another and at the indicated conditions. 

Let us now center our attention on the situation depicted in Figure 1. Init C, and 
Cp are cylinders fitted with weightless and frictionless pistons. The connection 
of the cylinders to the reaction box Cp includes valves V, and Vg —initially 
closed— as well as semi-permeable membranes Mg and Mg; the former 
permeable only to ideal gas a; the latter only to ideal gas f. The reaction box Cp 


ay <a 


contains a vast amount of a mixture of gases a and f at their respective 
equilibrium pressures P, og, Pp eq, at temperature 7, i.e., in Cp: 


a(Py,eq T)= B(Pe.eq T) (45) 

The pressure inside the box is constant and equal to Peg = Py,eq + Peeq. The 
pistons have been properly coupled to a work or mechanical reservoir (mr) -not 
shown- whose function is to serve as deposit or provider of the work either 
produced or demanded by the operation. Initially, while cylinder « contains 1 
mole of a at 1 bar and temperature 7, cylinder f is empty, with its piston in its 
extreme leftmost position. At this point the pressures keeping these pistons in 
place are 1 bar, and Pp og, respectively. The equilibrium box is in thermal contact 
with a heat bath or heat reservoir (bath) —not shown- of temperature 7: 


Gases a. and f in the cylinders, as well as their equilibrium mixture in the reaction 
box, constitute the thermodynamic system (sys); with the heat and mechanical 
reservoirs constituting the surroundings. Combined, system and surroundings 
constitute the universe of the process. The mechanical device constituted by 
reaction box, cylinders, pistons, valves, and semi-permeable membranes is 
called the van ’t Hoff Equilibrium Box (EB). 
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Figure 1. The van’t Hoff Equilibrium Box. 


The form of the first law to be here used to describe and quantify the energy 
exchanges of the system with its surroundings, is the following: AU = Q — W. 
In it, heat absorbed, and work done by the system are positive magnitudes; their 
inverses, negative. It is understood that in the universe just described, the only 
bodies capable of exchanging heat are the heat bath and the system, in which 
case heat lost by one, is heat gained by the other. Likewise, the only bodies 
capable of exchanging energy in the form of work are the system and the work 
or mechanical reservoir, in which case, work done by one is work done on the 
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other. It will be further assumed that the reader is knowledgeable about the 
thermodynamics of ideal gases. 


7.2. The four processes taking place in the van’t Hoff equilibrium box 
Process 1: The mole of gas a in cylinder Cy is expanded isothermally (7) and 
reversibly from 1 bar to Py eg. At the end of this process the pressure keeping 
this piston in place is P, 4. In accordance with the isothermal nature of this 
process, and the sole dependence of the internal energy of ideal gases with the 
temperature, we will have: 
AUsys4 = 0 (46) 

As indicated by the form adopted by the first law for this process, namely Qcy5 = 
W.ys, the reversible work transferred by the system to the work reservoir, comes 
from an equivalent amount of heat absorbed by it from the heat bath. It is with 
base on these considerations that the following equation will be written: 

—Qpatha = Qsys4 = Weys1 = —Wrra = —nhT In(P 2,eq) (47) 


In it, the last term on the right-hand side quantifies, as should be known, the 
reversible work of expansion produced by the ideal gas in its transit from P = 1 
to Pyeq. It should be noted that the argument of the logarithmic term in the 
previous expression is a quotient whose denominator, omitted here for reasons 
of economy of expression, is the unit valued initial pressure of the said 
expansion. 


At the conclusion of this process both valves are opened for Process 2 to occur. 


Process 2. Here n moles of gas a(Py,eq, 7) are introduced into the reaction box, 
and n moles of gas 6 (Pg eq,T ) simultaneously removed from it, 0 <n <$ 1. The 
respective representations of these two processes, given below: 


[Ne (e.aai T)] before adddition = [Ne (Piveai T)] after addition 
[ng (Pp.eq T) be fore extraction = [ng (Dp,eq) T)] after extraction 
allows us to realize that: 


Since the fugacities are the same in the cylinder as in the mixture, there 
is no change in free energy; the process occurs at equilibrium (Smith and 
Van Ness, pp. 409-410) 


The fact that no change in free energy means no exchange of energy in the form 
of work, leads us to: 


Wz =0 (49) 
The combination of the work associated to these two processes determines the 
work of the system for Process 2: 
Weys2 = Wa,2 a Wp = 0 (50) 
The sub-indexes here used refer to the changes experienced by the system, in 
the form of gases @ or f, in process 2 


That the temperature of each of these two gases remained constant along the 
processes just described, means that they experienced no change in their internal 
energies, i.e., : 

AU sy 52 = AUg2 + AUg,2 =0 (51) 


The substitution of the two previous results in the equation for the firs law, 
AUsy5 = Qsys — Weys, leads to the conclusion that: 


Qsys,2 =0 


The previous results evince the fact that no heat exchange takes here place 
between the system and the heat bath, that is: 


Qhath,2 = Qsys,2 =0 (52) 
Additional comments on this process can be found in Section 15. 
At the conclusion of this process the valves are closed. 


Process 3. Following the introduction/removal process just described, the 
following reaction takes place: 


na(PreqT) > MB(Pgeq:T) (53) 

The fundamental assumption behind the modus operandi of the van’ Hoff 
equilibrium box in general, and of reaction (53) in particular, or of the inverse 
of this reaction, is that even if the disequilibrium created on the vast amount of 
reaction in the reaction box by the addition of » moles of a, and the removal of 
n moles of 2 produces no measurable effect on the respective partial pressures 
or concentrations of these two chemical species, it will be however (this 
disequilibrium) the reason for the occurrence of reaction (53) through which 
equilibrium in the box will be restored. 


In Fermi’s words: 
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We assume that the content of the large box is so great that the relative 
changes in concenwation resulting from this inflow of gases is negligible. 
The concentrations of the gases during this process, therefore, remain 
practically constant ... Moreover, since the temperature and volume of 
the large box do not change, the chemical equilibrium of the gases in the 
box readjusts itself... (pp. 103, 105). 


The fact that in this reaction butane at its equilibrium conditions is transformed 
into isobutane at the same conditions is the reason why this reaction is said to be 
taking place ‘at equilibrium.’ The fact shown above in eq. (27), that at 
equilibrium the chemical force is zero, means that the occurrence of reaction 
(53) will take place without any exchange of work with its work reservoir: 


If fog = 0, then W,,,.3 = 0 (54) 


At the light of the previous result, the expression of the first law for reaction 
(53), that is for the system in process 3, namely AUsys3 = Qsys3 — Wsys,3> 
reduces to: 


AUsy53 = Qsys,3 (55) 


The previous equation shows that the whole of the internal energy or enthalpy 
change of reaction (53), which is the energy coming from the breaking and 
forming of chemical bonds defining this chemical process, appears as an 
equivalent amount of heat. 


The situation to be noted at this point in process 3 is that there exists no avenue 
open, no mechanism at play which could take this energy (the 
thermodynamically mandated fraction of if, that is), as it is being made available 
by reaction (53), and transform it into work. The only path open to it is the one 
leading to the heat bath, and it is this fact what defines the way through which 
the said fraction of this energy is to be converted into work. 


Let us acknowledge now the following two facts: 1) The internal energy, and 
the enthalpy, of an ideal gas are sole functions of its temperature (Denbigh, p. 
113). The previous statement implies that said functions are independent from 
the pressure; 2) Any ideal-gas phase, isothermal, and isobaric reaction taking 
place with constant total mole number, An = 0, like the one being here 
considered, is also isochoric; and on reason of this, for any such reaction the 
internal energy and enthalpy changes become identical: AU = AH — PAV = 
AH — RTAn= AH. 


It is with the previous notions in mind that the following equation will be written: 
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AU sys,3 = AH sys,3 > AH[na(Peeq T) > np (Ppeq T)| = n4H[a(1bar, T) 
— B(1 bar, T)] = ndH?® = nAU? (56) 
From which the following equation will be excerpted: 
AUsys,3 = nAU® = nAH° (57) 


The combination of equations (55) and (57) leads us to the following equation 
in which we have included the change experienced by the heat bath: 


n4H° = Qsys,3 = —Qhath3 (58) 
Process 4. Here, the amounts of a and £ in Cy and Ce at the end of process 3, 
namely (1—n)a@ (Poeq,T) and nB (Pg.eqT), will be now compressed 
isothermally (7) and reversibly to a final pressure of 1 bar. The amounts of work 
involved are the following: 
Waa a Qe,4 = —Qbath,a,4 = (1—n)RT In Pa,eq (59) 
Wea = Q¢,4 = —Qbath,p.4 = NRT In Py eq (60) 
Combined, the two previous equations quantify the changes experienced by the 
system and the bath in this process in the following manner: 
Wsys,4 = Was + Wea = Qsys,4 = —Qpath 
= (1 —n)RT In Peg + MRT In Pp og (61) 
The isothermal nature of these two processes means that they occur without 
changes in their respective internal energies. This fact is implicit in the equality 
between work and heat evinced by the previous three equations. Therefore: 
AU sys.4 = 0 (62) 
7.3 A summary of changes 
7.3.1 The reaction Box 


No change remains in the reaction box once the previous processes have 
concluded. The original changes consisting in the addition of m moles of a and 
the simultaneous removal of n moles of 8 were offset by the occurrence of 
reaction na(Paeq T) — NB (Ps .eq,T) through which the n moles of a inserted 
were consumed in order to replenish the 1 moles of £ extracted; with the 
constancy of its temperature assured by the action of the heat reservoir. 

7.3.2 The Cylinders 


Cylinders a and £ do experience changes in their respective conditions because 
of these four processes. Thus, while the former originally contained 1 mole of 
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a(1bar, T); at the end it contains only (1 — n) moles of this gas at these same 
conditions. The latter, originally empty, is found at the end containing n moles 
of 6(1bar,T). Note then that the previous concatenation of processes has 
managed to reversibly transform the original mole of a(1bar, T) into pure a and 
pure £ in the amounts and conditions given by (1 —7n)a (1bar,T), and 
np (1 bar,T). The former occupying cylinder a and the latter cylinder £. In 
other words, these four processes combined to bring forward the reversible 
standard reaction given by equation (44). 
7.3.3 The System 
The work performed by the system along these four processes, quantified in 
equations (47), (50), (54), and (61), amounts to: ; 

Woys = —RT I Peg +0 +0+(1—N)RT IN Pyeqg + MRT IN Pp oq (63) 


Upon simplification, this equation reduces to: 
P 
Weys = MRT In (2) (64) 
B a,eq 
Recognition of the fact that the pressure quotient in the previous equation is the 


equilibrium constant K for reaction (1) at temperature T, combined with the fact 
that: 


AG° = —RT InK (65) 
Allows writing equation (64) as follows: 
Woys = —nAG° (66) 
The heat exchanged by the system in each of these processes has been evaluated 
in equations (47), (52), (58), and (61). It amounts to: 

Qsys = —RT In Pajeg +0 + NAH? + (1 ~n)RT In Pyeg + NRT ln Pgeq (67) 
Simplification of this equation, followed by its combination with equation (65) 
produces: 

Qsys =n(4H°-AG°) (68) 
The defining equation of the Gibbs energy, G = H —TS, allows writing the 
previous expression as 
Qsys =NTAS 2 (69) 


The internal energy changes experienced by the system are given in equations 
(46), (51), (57), and (62). Combined, they amount to: 
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AUeys =0+0+nAH°+0=nAH? (70) 


The final form of the previous equation comes through the substitution in the . 
equation for the first of the results shown in equations (66) and (67), i.e., 


AU; =nAH® = Qsys — Weys = n(TAS° + nAG°) (71) 
7.3.4 The mechanical reservoir 


That in the universe of the equilibrium box the only bodies capable of 
exchanging energy in the form of work are the system and the work reservoir 
(mr) means that work done by one is work done on the other, and vice-versa. In 
other words, W,,; = —Wy,. According to this, at the end of the four processes 
above described, the work in the work reservoir finds, in attention to eq. (66), 
the following quantification: 


Weys = —Wrr = —nAG° (72) 


It should be kept in mind that while AG° in the previous equation is a constant 
which stands for the difference between the standard Gibbs energies of 
formation of 1 mole of isobutane and 1 mole butane, as they appear in reaction 
(1), both work terms there appearing W.y;, and Win,, refer to the conversion of 
n moles of butane into isobutane, 0 << 1. As expressed in the previous 
equation, both these terms are functions of 7. 


7.3.5 The heat reservoir 


That in the universe of the equilibrium box the only bodies capable of heat 
exchange are the system and the heat bath means that heat gained by one is heat 
lost by the other. According to this, and in attention to eq. (69), we can quantify 
the heat exchanged by the heat bath with the following expression: 


Qhath = —NTAS° — (73) 


All the results obtained in the previous discussion agree with thermodynamic 
theory. 


Regarding the standard isomerization reaction of butane into isobutane, we have 
this exothermic reaction making available in process 3, an amount of heat equal 
to n4H° = n+ 8368.4 J which is originally transferred to the heat bath. A 
portion of this energy in the amount of n4G° = n+ 3747.1J leaves the heat bath 
to be transformed into an equivalent amount of work that ends up in the work 
reservoir. The difference, in the amount of n+ TAS® = n(298.15- 15.5) =n- 
4621.3 J remains, as unavailable energy, in the heat bath. 


-72 « 


8. Back to the beginning 


It was through the concatenation of the four processes just described that the 
following reaction, as represented in eq. (44), took place: 
1a (g,1 bar,T) > [(1 — n)a@ (g,1 bar, T)] + [nB (g, 1bar, T)] 

This new concatenation of the four processes to be described below, will be 
responsible on the other hand, for the inverse change, namely: 

1a (g,1 bar,T) — [(1 —n)a@ (g,1 bar, T)] + [nB (g, 1bar, T)] (74) 
The processes and associated changes subsumed by it, will be identified with 
Roman numerals. When necessary, these processes/changes will also bear the 
subscript 6 > @ to distinguish them from those associated to the previous 
concatenation, which, when needed, will be identified by the subscript a > B. 
The new concatenation starts at the conditions we find the universe (system and 
surroundings) at the conclusion of the previous one. Here, the achievement of 


the intended change represented in eq. (74) will include, among other processes, 
the occurrence of the inverse of reaction (53), that is, of reaction: 


nB(PpeqT) > na(PxeqsT) (75) 
The description of these four new processes is given below. 
Process I: Gases a and & in the amounts and conditions they are respectively 
found in cylinders Cy and Cg at the conclusion of the previous concatenation, 
namely (1—n)a(1bar,T), and nf(1bar,T), will be isothermally and 
reversibly expanded to respective pressures Pyeq and Pgeg. The following 
changes accompany these expansions: 


AU sys,1 =0 (76) 


Wsys,a,I = Qsys,a,I = —Qhath,a,1 = —Warat 
= —(1—n)RT In(Pyeq) (77) 
Weys.p,t = Qsys,6,1 = —Qbatn,g,1 = —Wmr,p1-= —NRT In( Pp eq) (78) 


The fact that Woy.67 = Woys.a,1 + Weys,g,1; leads us to the following expression: 


Weys,r = Qsys,1 = —Qpatht = —Wrrt 
=—-(1- n)RT In(P2,eq) — nRT In(Pp,eq) (79) 


After rearrangement and simplification, the previous equation takes the 
following form: 
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Wsys,t = Qsys,t = —Qhati,t = —-Wnrt 
= —RT In(Pyeq) + MRT In (Po,eq) — NRT In (Pe,eq) 


= Foeq 

= —RT In (Py¢q) + NRT In (——) 
Pp,eq 

= —RT In (Pyeq) — NAG _,e (80) 


Process II. The fact that the addition/removal process taking place here is 
essentially the same as that discussed in process 2, allows us to identify the 
introduction of n moles of gas 6 (Pp eq, 7) into the reaction box, and the removal 
from it of 7 moles of gas @(P,,eq,T) from it, as processes taking place with no 
change in Gibbs energy, and if so, without any exchange of energy in the form 
of work. Therefore: 

Wan = 9 (81) 

Wen =0 (82) 
Note here that at the conclusion of this process we will find that while C, 
contains 1 mole of a, Cg is, on its part, empty. 


The combination of the two previous equations produces: 

Weys tl = 0 (83) 
From the isothermal nature of these two processes, we get: 

AU gys,11 =0 (84) 
And from the first law, in the form Qsys.17 = AUsys,11 + Weys,1, the following 
result is produced: 

Qsys1 =0+0=0 (85) 

Process III. It is in this process where, via eq. (75), the 7 extra moles of f in the 
reaction mixture are consumed to offset its 7 moles of a deficit. The fact that it 
is a reaction taking place ‘at equilibrium,’ means, as already noted regarding 
these reactions, that no work is exchanged by it, if so, then: 

Weys tll = 0 (86) 
In terms of the previous result, the first law reduces to: 


AU sys 111 ae Qsys,HI (87) 
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The fact made evident by a simple inspection of eq. (56) that AU sys 117 = 
—AU sy 3 allows us, on reason of eq. (57), to re-express eq. (87) in the following 
manner: 


AU sys ii = Qsys,tII = AH3 +a (88) 
The previous equation shows that along its occurrence, reaction (75) takes away 
from the heat bath an amount of heat equal to AH, Bae identical in magnitude to 
that transferred to the heat bath by reaction (53) in process 3. 
The following quote from Clausius serves to set processes 3, and JII, in the 
proper thermodynamic frame: 
... let us consider a cyclical process, which is such that a body passes 
through a series of changes of condition and at last returns to its original 
state. This variable body, if placed in connection with the heat reservoir 
to receive or give heat, must have the same temperature as the reservoir; 
for it is only in this case that the heat can pass as readily from the 
reservoir to the body as in the reverse direction, and if the process is 
reversible it is requisite that this should be the case. (p. 106) 
Process IV. At this point the 1 mole of @ in cylinder Cg is compressed 
isothermally and reversibly from Py,¢g to 1 bar. The changes associated to this 
process are the following: 


AU sys,1v = 0 (89) 
Wsys.iv = Qsys,1v = —Qpatniv = —Wmriv = RT In Poyeq (90) 
8.1 Summary of Changes for the inverse path. 


8.1.1 Internal energy: From equations (76), (84), (88), and (89) we can write, 
for the internal energy change of the system along the inverse process (inv), the 
following equation: 


AU sysinv = 0 +04 NAHB gq + 0 = NAHB. (91) 
8.1.2 Work: The work exchanged by the system is determined by the 
combination of equations (80), (83), (86), and (90): 
Wsys,inv = —RT In(e,eq) — nNAGPB +a + RT In@aeq) = —n AGB sa (92) 


8.1.3 Heat: The heat exchanged by the system with the heat bath can be obtained 
from equations (80), (85), (88), and (90): 


N 
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Qsys.inv == —RT In(Pa,eq) = NAGR +a + 0 + NAHB +a + RT In(De,eq) 
= n(AHp +a - AGB 4a) = nTASB 4a (93) 
9. Summary of results for the cycle defined by the combination of the two 
previous concatenations 


9.1 Internal energy: From the combination of equations (70) and (91) we get: 
AUsyscycie = AHg-4p + AHR = 9 (94) 


The equality to zero of the previous equation stems from the state function nature 
of the enthalpy, which asserts that AH? _,g = —4Hp 4a: 


9.2 Work: The work exchange by the system and the work reservoir, calculated 
with equations (72) and (92), leads us to the following result: 


Wsys,cycle = —Wnrcycle =—nd GBoa re nAGy.p = 0 (95) 
The zero in the previous equation finds the same explanation as that given for 
eq. (94), 
9.3 Heat: The heat exchanged by the system with the heat bath, calculated with 
equations (69) and (93): 


Qsys,cycle = —Qhathcycle = NT AS¢ 4g + NT ASB +a =0 (96) 
Note that here, again, nTAS¢_,2 = —nTASB_.q. 


Equations (94), (95), and (96), make evident the fact that at the completion of 
the cycle, no change remains in the system, the bath, or the work reservoir. At 
the end of this cycle, the universe has returned to its original condition, and no 
evidence remains of any change having ever taken place; which is the essential 
meaning of the reversibility concept, at least in the way it was originally 
understood by Carnot (pp. 11, 19). Said restoration implies the reversibility of 
both, the forward, and reverse processes. 


The introduction into the reaction box, in the first concatenation, of only a 
fraction n of the one mole of available butane obeyed to the fact that the results 
thus obtained could be concatenated with the mixing process to be described 
below, in such a way that the result of this combination could be compared with 
that of the chemical-potential centered analysis represented by eq. (43) 


10. The mixing process 


It was Lord Rayleigh the first to recognize that work might be obtained by the 
isothermal mixing of gases. In a paper titled “On the work that may be gained 
during the mixing of gases,” published in 1875, he demonstrated, in the words 
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of Darrigol (2018), that “... the maximal work that can be obtained by mixing 
(isothermally) two gases initially occupying the separate volumes V, and V, and 
finally sharing the volume V,+V2 is equal to the work produced by the 
expansion of the first gas alone from V; to V; + V2 pus the work produced by the 
expansion of the second gas alone from V2 to V; + V2. The maximal work is 
reached when the mixing is reversible.” 


The following discussion pertains the isothermal and reversible mixing of the 
product of the first concatenation, that is, of the (1 — n)a@(1bar,T) moles of o 
in cylinder Cy, and the nf (1bar, T) moles of f in cylinder Cg. Apart from saying 
that this isothermal, reversible, and work producing mixing process takes place 
in a device in thermal contact with the same bath, and properly coupled to the 
same work reservoir used in the van’ Hoff operation; we will add that excellent 
discussions about this issue are available in the thermodynamic literature (Klotz 
and Rosenberg, pp. 271-273; Schmidt, pp. 135-138). In them we find the mixing 
work given by the following equation: 
Weysmix = —(N@RT In Xq + ngRT In xg) (97) 
In terms of the notation which based on the degree of advancement €, was 
introduced in eq. (17) of Section 2, making ng = xg = pg = §, and ng = X%q = 
Pe =1-—€, the previous equation becomes: 
Woysmix = —RT[(1 — €) In(1 — &) +€ In &] (98) 
Applying the same change of variable to eq. (66), which gives the work 
produced by non-mixing standard reaction (44) in the first concatenation, we 
get: 
Weys = —fA Gop (99) 
At the end of this mixing process, we find that the total work produced by the 
system amounts to: 


Wsys,tot = —$AG3_,g — RT[(1 — §) In(1 — §) + § Ing] (100) 
In terms of free energy, the two previous equations find, on reason of eq. (5), the 
following equivalent expressions: 
AGgys = EAGE_,g (101) 
AGsystot = FAGZ_,p + RT[(1 — §) In(1 — §) + F Ing] (102) 


A simple inspection will reveal the identity existing between equations (43) and 
(102), confirming this way the equivalence of the usual thermodynamic 
procedure with that of the van’t Hoff equilibrium box. 
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11. Discussion 


11.1 The previously noted fact that reaction (44) as it takes place in the van’t 
Hoff equilibrium box transits from an initial state represented by 1 mole of pure 
butane in cylinder C, to a final state represented by (1 — n) moles of pure butane 
in the same cylinder, and n moles of pure isobutane in cylinder Cg, both at 1 bar 
of pressure at the selected temperature, is what makes of this reaction a standard 
reaction: 


The word ‘standard refers not to any particular temperature, but to unit 
pressure of 1 bar for each of the pure reactants and products.” (Denbigh, 
p. 148) 


The ‘pure’ condition means, of course, unmixed, 
The standard path for reaction (44), a non-mixing path, is thus the one defined 
by eq. (101). The fact indicated by Table 2 that AGZ_,, = —3747.1J mol™’, 
allows to write this equation in the following manner 

AGsys = —3747.1n = —3747.1€ (103) 


Any point on the negatively sloped line defined by this equation on coordinates 
AGsys — vs.—§ corresponds this way to a different run of reaction (44); Each 
of these runs, characterized by a given n or &, producing a different change in 
the Gibbs energy of the system, and a different work output. The fact that for all 
n different from zero, 0 < n < 1, it is true that AG,,,, < 0, means that reaction 
(44) is spontaneous for all these n. 


The fact that the linear path of reaction (44) is negatively sloped or, equivalently, 
a Gibbs energy decreasing path, means that any number of moles of a@ introduced 
into the van’t Hoff box will be totally converted into the stoichiometrically 
corresponding moles of /. 


It is interesting to note here that this devoid of mixing, and full-conversion 
standard reaction comes out of the equilibrium mixture of the same reaction, a 
reaction which precisely on reason of having taken place along a mixing path, 
was incapable of complete conversion. 


11.2 As proved above, the mixing of the products of reaction (44) transforms the 
path of reaction (44) from one described by eq. (101) to one described by eq. 
(102). The graph of this equation takes the familiar form of the concave up curve 
with its minimum corresponding to the equilibrium mixture. 
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Of the two sources capable of providing energy for reversible reaction (1) to 
transform into mechanical work, it is the one coming from the chemical process 
itself, that is, the one associated to the breaking and forming of bonds (nAG@.,2) 
the one with which the standard reaction deals with. The ‘actual’ path of a 
chemical reaction, deals with both, the previous one, and the one originating in 
the mixing process (Ng RT In ng +7gRT In ng). 


11.3 In thermodynamics, reversible is any process driven by an internal force 
infinitesimally larger than an external resisting force of the same nature. Thus, 
the reversible expansion of a gas requires the pressure of the gas being 
infinitesimally larger than an external, resisting pressure; For example, 
maximum work is obtained from a chemical reaction taking place in an 
electrochemical cell working “against an external voltage lower than its own by 
no more than a differential amount” (Pimentel & Spratley, p. 112). 


The fact that in chemical reactions with non-ionic mechanisms, say the gas- 
phase reaction above considered, it is apparently impossible to oppose a force 
of the same nature to the one driving the reaction, that is, to the one which called 
the ‘chemical force’ and represented by f, finds definition in terms of the 
difference of the stoichiometrically weighted chemical potentials of reactants 
minus products, seems to be the reason why the only practical path open to the 
reversibility of these reactions is the one represented by the van’t Hoff 
equilibrium box in which the reaction of interest is seen taking place ‘at 
equilibrium.’ The reason why this is so simple: A reaction taking place at 
equilibrium does so with a chemical force equal to zero. In this case the opposing 
force is redundant, there is just no need for it given that there is nothing for it to 
oppose. 


11.4 We have already proved the usefulness of equation 
dW = fdé for the quantification of the work produced by a reversible reaction 
along any given transit ¢, — €,. The interesting thing to note in this regard is 
that the work coming out of the van’t Hoff equilibrium box does so without the 
intervention of this force. It could not have been otherwise given that the reaction 
in it is taking place ‘at equilibrium,’ condition where this force becomes equal 
to zero. What we find instead is that the work coming out is actually PdV work. 
The ‘interesting’ qualifier used above pertains the fact that the work coming out 
of a reaction taking place without any change in volume, such as reaction (1), 
finds evaluation via dW = PdV 
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11.5 The van’ Hoff equilibrium box is a tool devised for the theoretical analysis 
of chemical reactions, with no possibility of practical realization. Even if you 
succeed in creating infinitesimal gradients, you will have to wait an infinite 
amount of time for the process to conclude. If these two restrictions were not 
enough, the case of the equilibrium box faces the additional problem of 
availability of the required semi-permeable membranes. These are the reasons 
why instead of using the equilibrium box to produce, for example, mechanical 
work out of the Gibbs energy change associated to the reversible oxidation of 
methane in a van’t Hoff engine, we are forced to use its enthalpy change for this 
purpose, and this is done by burning it and feeding the heat thus produced to a 
less efficient, finite-time, thermodynamic cycle (Denbigh, pp. 72-76), a path 
associated to an undesirable subproduct: pollution. Let me note in this regard 
that the theoretical foundation provided by a corrected and extended version of 
quantum chemistry called Hadronic Chemistry has made possible the 
development of new clean energies and fuels (Santilli, 2001). This discipline, 
composed of three branches: Isochemistry, directed to the representation of 
closed isolated systems such as molecular structures; Genochemistry, directed 
to the description of “... irreversible chemical processes with unrestricted 
interactions, as expected in chemical reactions.” (ibid, xl); and Hyperchemistry, 
directed to the representation of “... irreversible biological structures with multi- 
valued unrestricted internal processes.” (ibid, xli), has developed new models 
for the hydrogen and water molecules in agreement with “... conceptual, 
theoretical, and experimental evidence.” (Santilli & Shillady,1999; Santilli & 
Shillady, 2000, p. 81) 


12. The van’t Hoff equilibrium box as the chemical counterpart to Carnot’s 
heat engine 


Carnot’s reversible heat engine produces work out of heat flowing from a hot 
reservoir of temperature T;, to a cold reservoir of temperature T,, and does it by 
introducing between these reservoirs what is called a Carnot engine, basically a 
working substance, say one mole of an ideal gas, which contained in a cylinder 
fitted with a frictionless and weightless piston properly connected to a 
mechanical reservoir, is subject to the cyclical concatenation of four reversible 
processes: 1) an isothermal and reversible expansion at the temperature of the 
hot reservoir. The conditions of the gas in this process change from (T}, V4) to 
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(Tp, Vg); 2) an adiabatic and reversible expansion lowering the temperature of 
the gas to that of the cold reservoir; the conditions of the gas transit here from 
(Tr,Vg) to (Tc,Vc); 3) an isothermal and reversible compression at the 
temperature of cold reservoir, changing the conditions of the gas from (T;, V,) 
to (T,;,Vp); and finally, an adiabatic and reversible compression which starting 
at (T;,Vp) ends up with the gas at the same conditions it originally had at the 
start of the cycle, that is (T;,,V,). The previous four processes, in the order 
mentioned, will also, for reasons of economy of expression, be designated as 
AB, BC, CD, and DA. Each of the engine’s work producing cycles starts with 
the hot reservoir transferring to the gas along process AB, an amount of heat Qp. 
The cyclical work output will, on its part, be designated as W. 


With the previous antecedents at hand, let us recognize the following salient 
operational characteristics of this engine: 


1. Combined, adiabatic and reversible processes BC and DA make no 
contribution to the work output of the cycle. The reason being that the work 
produced by the former is exactly the amount consumed in the occurrence of the 
latter. For adiabatic processes, characterized by Q = 0, the first law reduces to 
AU = w. For ideal gases, this equation can be written as AU = C,,AT. Applied 
to processes BC and DA, the previous equation takes the following forms: 
AUge = Cy(Te — Th) and AUp, = C,(T, — T;). These equations, combined, 
lead to AUge + AUp, = 0, and on reason of the just noted identity between 
internal energy change and work, also to Wgc+Wp, =0. Even if not 
contributing to the work output of the cycle, isentropic processes BC and DA are 
the ones making possible the cycle. 


2) The previous result leads to the realization that the work output WW of the cycle 
is the sole matter of processes AB, and CD. In the former, an isothermal and 
reversible expansion of the ideal gas at the temperature T;, of the hot reservoir, 
the whole of the heat Q;, received by the gas from the hot reservoir is transformed 
into an equivalent amount of work w; which ends up in the work reservoir. Once 
the temperature of the gas has been taken from T;, to T, by process BC, the gas 
is subject to an isothermal and reversible compression at the temperature of the 
cold reservoir, T,. The work consumed to produced this change, w,, represents 
a portion of the one previously deposited in the work reservoir by process AB. 
In line with the mandate of the first law, as it applies to this process, the work 
spent on it ends up as an equivalent amount of heat Q, in the cold reservoir. 
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3. The difference of the temperatures at which expansion AB and compression 
CD take place, explains, among other things, the partial conversion into work of 
the heat released by the hot reservoir. 


For isothermal and reversible processes AB and CD, the following equations can 
be respectively written: wy, = Qn, =RT;,In(Vg/V4) and w.=Q,= 
RT-In(Ve/Vp) , in which case the work output of the cycle will be given by 
W = wy — We. If we now bring here the equation relating the initial and final 
states of isentropic processes BC and DA, namely Vg/V, = Ve/Vp, and with it 


write Wo = Q, = RT, In ve , and once this is done, we take the quotient w;,/W,, 
A 


we will get w,/W, = T,/Te. This last expression, on substitution of w, = w;, — 
W, becomes (w;, — W)/w;, = T-/Th, which on further algebraic rearrangement 
produces: 


Ww w heat 


Wh = On = Ncarnot = a ae 

The said difference in temperatures, even if restricting the full conversion of heat 
input into work output, is the sine qua non condition that opens the door to the 
continuous transformation of heat into work in a heat engine. When the heat 
engine returns to its initial condition to start a new cycle, the work reservoir is 
found holding an additional amount of work W = YegrnorQp to the one it holds 
at the conclusion of the previous cycle. Let us now note that an isothermal and 
reversible ideal gas expansion has the capability to transform its heat input into 
an equivalent amount of work. This, however is a one-time event. Any attempt 
to construct a work producing cycle by concatenating this expansion with its 
inverse will fail on reason of the fact that at the start of the new cycle no work 
will be found in the mechanical reservoir, the reason being that in the isothermal 
and reversible compression needed to produce said cycle, the whole of the work 
initially outputted will be consumed making it a zero efficient cycle. No work 
can come out of an isothermal Carnot’s cycle. 


13. The van’t Hoff chemical-potential engine in comparison to that of 
Carnot 

The following arguments will refer to the operation of the van’t Hoff equilibrium 
box, when instead of the introduction of n@ (pg,eq, T) moles of butane to the vast 
amount of equilibrium mixture in the reaction box, and the simultaneous 
removal from it of n6(pgeq,T), where 0 <n <1, we introduce to the said 
equilibrium mixture 1 mol of butane, and remove from it 1 mol of isobutane, 
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both at the previously indicated conditions. The benefit of this procedure is the 
simplicity of description afforded with it. 


The form adopted by the four processes constituting the operation of the van’t 
Hoff engine in this new case are the following: 


Process i: The isothermal and reversible expansion of the 1 mole of butane in 
cylinder Cy from 1 bar to pg,eq. From equations (46), and (47), this last one via 
the substitution n = 1, we get: 


AUsys,i = 0 (104) 
—Qbatni = Qsys.i = Weys.i = —Weysi = —RT In(Po,eq) (105) 
Process ii: The introduction of 1 mole of @(Pgeq,T) to the vast amount of 
equilibrium mixture in the reaction box, and the simultaneous removal from it 


of 1 mole of 6 (Pg,eq T)i For the reasons already given in regard to processes 2, 
and II, we will have that: 


Wa2i = 0 (106) 
We,2i = 9 (107) 
Weys.2i = Wa,2i + Weri = 0 (108) 
AUsys2i = AUg,2i + AUg.2i = 0 (109) 
Qpath,2i = Asys,2i = 0 (110) 
Note that at the end of this process while cylinder Cy is empty, Cg contains 1 
mole of B (pg eq, T) 


Process 3i: Equilibrium is restored in the vast amount of reaction inside the 
reaction box via the occurrence at equilibrium of the following reaction: 


1@(Pe,eq T+) 7 1p (Pp,eq T) 


The ‘at equilibrium’ nature of the previous reaction means, as already discussed, 
that: 


Weys,3i = 0 (111) 
From equations (57) and (58), via the same substitution, we get: 
AUsys,34 = 4U° = AH® (112) 
AH® = Qsys,3i = —Qpath,3i (113) 
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Process 4i: The mole of £ in cylinder Cg is compressed isothermally and 
reversible from pp,eq to 1 bar. With base on (60) and (62) and the previously 
noted substitution, we get: 


Wsys,ai = Qsysai = —Qhatnsi = RT In Ppeq (114) 
AUsysai = 9 (115) 


14. Van’t Hoff’s chemical-potential engine compared to Carnot’s heat 
engine 

With the two previous sections in place, the following statements will be 
advanced: 


1. The obvious difference between Carnot’s and van’t Hoff’s engines is that 
while the former exploits the work producing potential subsumed by the transfer 
of heat from a higher to a lower temperature, the latter produces work out of the 
“_.. transfer of chemical matter between two different chemical states described 
by chemical potentials.” (Gill, p. 507). The fact that the purpose of one and the 
other is the same, namely the production of work, is the reason why the van’t 
Hoff equilibrium box is being here referred to as an ‘engine.’ 


Just like no work is outputted in the absence of a temperature gradient in 
Carnot’s engine, neither is, in reference to a chemical reaction, when the value 
coming out of the combination of the chemical potentials of the reactants, 
properly weighted by their respective stoichiometric coefficients turns out to be 
identical to that of a similar combination of the chemical potentials of the 
products. In both cases the reason is the same: from a system in equilibrium 
(thermal in the former, and chemical in the latter) no work production is 
possible. 


2. A closer look at equations (105), (108), (111), and (114) allows us to realize 
that the work output of van’t Hoff’s engine finds quantification in precisely the 
same way it happens in Carnot’s engine, that is, in terms of the combination of 
the work produced by an isothermal and reversible expansion, and the work 
consumed by an isothermal and reversible compression. The difference being 
that while in Carnot’s engine they take place at different temperatures, in van’t 
Hoff’s engine they take place at the same temperature. 


3. In a heat engine it is the no-work-contributors adiabatic and reversible 
processes the ones making the reversible cycle possible. They break the finite 
temperature gradient into an infinity of steps differing infinitesimally in 
temperature, and in doing so providing, via the working substance, a reversible 
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connection between them. Once a cycle concludes via process DA, the engine is 
ready to start a new one in which, again, one portion of its heat input Q;, will be 
transformed into mechanical work, with the rest discarded into the cold reservoir 
as the amount of heat Q, = Q, — W, and so on for each other cycle. 
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Figure 4. The reversible transit of 1@(1 bar, T) — 18(1 bar, T) in a van’t Hoff 
engine is here shown taking place via the concatenation of work-producing 
isothermal and reversible expansion 1, work-less processes 2 and 3 jointly 
representing the ‘at equilibrium reaction’, and work-consuming isothermal and 
reversible compression 4. 


In van’t Hoff’s engine, on the other hand, it is the ‘at equilibrium reaction’ 
—which for the sake of this argument will be taken as the combination of 
processes 2i and 33 (after all, reaction 3i won’t take place without the previous 
occurrence of process 27) the one which: 1) in combination with work-producing 
process li, and work consuming process 4i, provide the reversible path to the 
transit of chemical matter from y2 to MB; and 2) by itself makes possible the 
cyclical operation of van’t Hoff’s engine. The reason behind this last assertion 
comes from the fact that the infinitesimal displacement which in the equilibrium 
condition of the vast amount of reaction in the reaction box is produced by the 
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addition to it of 1 mole of a, and the simultaneous removal from it of 1 mole of 
B, finds itself cancelled by the opposite-direction displacement through which 
said reaction recuperates its initial equilibrium condition via the consumption of 
the mole in excess of a in order to replenish the lacking mole of f. At this point, 
it should be noted, with equilibrium restored, the reaction mixture is in the 
condition to perform a new cycle which starts, as before done, with the addition 
of anew mole of « and the removal of 1 mole of f. Each cycle responsible for 
the output of an amount of work W = —AG@.z. 


Figure 4 provides a graphical representation of this operation. 


From equations (105) and (114), in combination with the following data, R = 
8.314] K"tmol™*, T = 298.15 K, Peeq = 0.18 bar, and pg eq = 0.18 bar, 
these last two items from Table 3, we get: 


Weysi = 1+ 8.314 + 298.15 « In(1/0.18) = 4250.67 J 
Wsysai = 1° 8.314 + 298.15 + In(0.82) = —491.92 J 


The algebraic summation of these two results produces the work output W of 
van’t Hoff’s engine: 


W = Weysi + Weys,ai = 4250.67 — 491.92 = 3758.5} (116) 


The number obtained above for W differs from the expected result, namely 
—AG° = 3747.1], in 11.4 J, or 0.3%, a reasonable difference explainable in 
terms of rounding off error. 


In a first approximation we could use, save further pondering on this issue, the 
procedure used in Carnot’s engine to determine the efficiency of the present 
operation of van’t Hoff’s engine; Thus: 


5 SOE hee, 7445 
Novan't Hoff ~ Weysi ~ 4250.67. oy) 


As evinced by the figure, process 4i, the one following the equilibrium 
restoration of the reaction in the reaction box, is non-spontaneous, and as such 
work-consuming. The removal of 491.92 J out of the 4250.67 J available in the 
mechanical reservoir are required to be spent to carry on this compression. The 
operation of the van’t Hoff box will be now complete, with an amount of 
available work in the mechanical reservoir of 3758.5 J. 


15. There is one aspect of the previous graph that needs clarification, and it is 


the one related to the fact that while the numbers used to draw the paths of 
processes 1i, and 4i in Figure 4 are those which quoted in Table 3 refer to the 
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reversible, mixing-included-evolution of reaction (1); the ones that should be 
defining said paths are those coming out of the isothermal and reversible 
expansion (process 17) and compression (process 4) of pure ideal-gases a and 
f, respectively. To understand that the values of Table 3 are common to both 
these instances starts by recognizing that the act of fixing the value of & to be 
substituted in the equations for the chemical potentials of these substances 
(equations (30) and (31)) defines a particular mixture of given partial pressures 
of butane and isobutane. In this mixture no more change takes place. Asserting 
the contrary implies, contrary to assumption, that ¢ can change. It is its 
unchanging nature what makes of it indistinguishable from a simple, non- 
reacting mixture of the said gases at the said concentrations or partial pressures: 
Pa = (1—€) bar for butane, pg = € bar for isobutane. But in this simple 
mixture of ideal-gases we have, in accord with J. W. Gibbs (p. 157), as well as 
Dalton’s law, that each ideal-gas is as a vacuum to the other, or, equivalently, 
each gas behaves as if it were alone (pure) in the total volume of the mixture at 
the temperature of the reaction. Thus, the chemical potential of butane in that 
reaction mixture characterized by € = 0.3 in which, in accord with eq. (17), 0.7 
moles of butane exert a partial pressure of 0.7 bar at the temperature of the 
reaction, is identical to the chemical potential of 0.7 moles of pure butane in the 
volume of the mixture at the same temperature. The total pressure exerted by 0.7 
moles of pure butane at temperature T when confined in a volume JV, of 
magnitude P, = 0.7 (RT/V), turns out to be identical to the partial pressure 
exerted by the same amount of gas in the reaction mixture at the same volume 
and the same temperature: p, = 0.7 (RT /V): 


Mg (pure, P, = 0.7 bar, T) = Mg (mixed, pg = 0.7 bar, T) 
= 2 + RT In(0.7) (118) 


In the words of Gibbs (p. 158) “... the potential 4, has the same value in the gas 
mixture and in the gas G, existing separately as supposed.” An exemplification 
of this notion can be found Castellan (p. 224-224). 


The argument leading to eq. (118) is the theoretical foundation of Smith and Van 
Ness explanation made in regard to the zero Gibbs energy change of process 2: 
If Ua(pure, Py,eq T) = Ha (mixed, PoeqsT), 


then Allsyso = Ua (mixed, Prot) — Ma(pure, Pusat) 
=0 (119) 
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Abstract 


There are two stages in the process of low-temperature nuclear fusion. 
The first stage corresponds to the bringing together of nuclei due to attraction to 
the negative charge in the internuclear space. This charge is formed by massive 
electron pairs located in a circular orbit. According to hadronic mechanics, the 
attraction of electrons (with opposite spins) is ensured by contact interaction. The 
second stage corresponds to the exchange of nuclei by pions, which begins when 
the nuclei approach each other to a distance equal to the strong interaction radius 
R,. To estimate the maximum value of R 5, the concepts of hadron mechanics 
about the x° - meson as a bound electron - positron (e" e") - pair are used. The main 
attention is paid to the visual interpretation of the electromagnetic interaction (€ 
e’) - pairs with nuclei, promoting the attraction of nuclei. 


Key words: low-temperature fusion of nuclei, contact interaction of electrons, 
strong interaction, pions, electron-positron pairs. 
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» 0. 
Introduction 


Currently, there is a lot of experimental evidence of low-temperature nuclear 
fusion, ranging from the simplest reactions of muon catalysis to the fusion of complex 
nuclei. As noted in [1], the synthesis of nuclei when a proton is captured by nuclei 
with charge numbers an orders of magnitude greater than unity, and, moreover, fusion 
reactions of complex nuclei at low temperatures should not occur due to tunneling. 
Such a synthesis requires bringing nuclei closer to distances at which the strong 
interaction of attraction between nuclei begins to exceed the Coulomb repulsion. Let 
us introduce the notation R; for the strong interaction radius. Obviously, the larger the 
value of Rs, the faster the convergence of nuclei required for the fusion reaction can 
be achieved, 


It is known that the strong interaction is short-range, i.e. Rs is orders of magnitude 
smaller than the atomic scale Ra~ 10°!’ m. Therefore, the “switching on” of the strong 
interaction must be preceded by the process of bringing nuclei together from Ra to 2d 
< Rs, where d is half the internuclear distance. According to the authors, this 
preliminary process is described relatively simply in the intermediate quasi- molecular 
state (IQS) model. The IQS diagram is shown in Fig. 1. 


ZAR. 


anh Ze 


2d 


> 


Fig.1. Scheme of the simplest model of a IQS ina state of electrostatic equilibrium 
(the R/d ratio is reduced compared to the calculated value V3, the CC activator 
contains two electron pairs) 


In this model, according to [1-4], the approach of nuclei is achieved due to the 
attraction of nuclei with a charge +q to a massive negative charge in the internuclear 
space. The negative charge is associated with a configuration in the form of compact 
massive electron pairs located in a circular orbit. This object, which catalyzes nuclear 
fusion, is called, for brevity, a CC activator. Pairs of electrons have opposite spins and 
are connected by an attractive contact interaction that dominates Coulomb repulsion. 
Nonpotential contact interaction is described by hadronic mechanics [5]. 


Estimation of the maximum value of the strong interaction radius 

We believe that the conditions for the first stage of nuclear convergence in the IQS 
model are met. Let us then show, following [1], that for the value of Rs we can use 
the estimate R; ~ 100 fm. 
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Let us recall that the mechanism of any type of interaction of objects that are 
not in contact with each other in quantum field theory is reduced to the exchange 
of virtual particles ~ the carriers of this interaction. In the case of strong 
internuclear interaction of interest to us, the mechanism first proposed by Yukawa 
[6] remains relevant. In this case, nucleons exchange virtual 2 - mesons, just as 
electric charges exchange virtual photons. 


Let us also recall that of the three x°, 2", 2 - mesons, the 2° - meson has the 
smallest rest mass (~ 135 MeV on the energy scale). Consequently, the estimate of 
the maximum value of the radius (Rs)max that interests us is asseciated with the 
exchange of nucleons by x° - mesons. The standard estimate of the interaction 
radius in quantum mechanics is given by 


R&c At=hc/2AE, (1) 


where c = 3 10° m/s and h =h/2 x= 1.05-10 ~™ Js and it is taken into account that 
the uncertainties in energy AE and time At are related by the Heisenberg relation 


AE At>h/2. (2) 


It is clear that at finite interaction radii, the maximum value of R corresponds 
to the minimum value AE of the order of the rest energy of the virtual particle, if it 
is neutral (x° - the meson is neutral). If virtual particles have an electric charge, 
then, in accordance with the law of conservation of electric charge, they should be 
born in pairs, so that the exchange of nucleons by x, x” - mesons is associated 
with a smaller radius R, compared to R, when exchanging 7° - mesons. 


Moreover, if we assume that, according to hadronic mechanics [5], the nm - 
meson is a bound state of an electron and a positron, then as the maximum value of 
R, we can formally take the value Rt ~ 10 ~° m. Indeed, the energy spectrum of 
the electron is positron pairs (associated with the x° - meson) lies in the range = (1 
- 135) MeV, that is, it is limited from below by the value of the positronium rest 
energy. Therefore, we take AE ~ 1.022 MeV as the minimum value. Then from (1) 
we obtain Ri ~ 100 fm = 10 ~ ' m. Note that the values of R% and the Compton 
electron wavelength are of the same order of magnitude. This indicates a 
characteristic intermediate scale (between the atomic and hadron scales). 


“Flickering mode” of the action of (ee’)-pairs on nuclei 


Let us recall that the classical model of positronium represents a pair of 
oppositely charged point particles (with masses equal to the rest mass of an 
electron) rotating in a circular orbit around a common center of inertia. In this case, 
the radius of the positronium orbit is approximately twice the Bohr radius, that is, 


«Os 


it has an atomic scale. Then a reasonable question arises: is it possible to compare 
the formal estimate of R} with an acceptable visual picture of the interaction of 
nuclei with the participation of positronium? Indeed, for nucleons located at 
distances d~R} the “exchange" by positronium, the size of which is of the order of 
R, >>R3, looks, to put it mildly, unconvincing. 


Nevertheless, since electromagnetic interaction dominates on the atomic 
scale, a completely natural and intuitive scenario can be proposed. Namely, we 
assume that under the conditions of the existence of !QS, the birth of electron- 
positron (e’e’) - pairs from the physical vacuum occurs predominantly with the 
localization of electrons in the central region of the CC activator. It is clear that 
such electron localization promotes the approach of nuclei not only due to their 
attraction to the electron, but also due to repulsion from the positron, whose 
orbit covers the region containing the nuclei. Of course, under the conditions of 
the existence of the IQS, the shape of the positron's orbit will differ from a circle. 
However, when viewed qualitatively, the shape of the orbit does not play a role. 


Note that for massive nuclei, the action of the electric field from a fast- 
moving positron can effectively manifest itself as a “flickering mode” that 
promotes the approach of one or the other nucleus, as illustrated by the diagram 
in Fig. 2. 
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Fig. 2. “Flickering mode” of the action of (e’e*)-pairs on nuclei. 
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Let's give an explanation. Rotation along a circular (or elliptical) orbit is 
uniquely associated with a pair of oscillations along orthogonal axes. The 
“flickering mode” can be interpreted as a consequence of the obvious axial 
symmetry under IQS conditions, which allows one to give preference to one of the 
vibrations, namely, vibration along the axis of symmetry passing through the 
nuclei. Indeed, in the presence of an axis of symmetry, the orbits of positrons that 
enclose nuclei and differ in the orientations of the orbital planes are equally 
probable. As a result, massive nuclei respond only to the axial component of the 
electric field of rapidly moving positrons, while the transverse component of the 
field (due to the rapid change in directions when the orbital plane rotates) is 
averaged. 


It is clear that at a high frequency of appearance and disappearance (e’e’) - 
pairs, the “flickering mode” becomes equivalent to the action of an electrostatic 
field, promoting the convergence of nuclei. In addition, as the process of 
“exchange” of nuclei (e’ e’) - pairs comes closer, more and more compact and more 
massive (e’ e*) - pairs are included up to the state of x° — meson. 


Conclusion 


The above considerations demonstrate how, already at distances of the order 
of 100 fm, interaction with the physical vacuum (PV) can occur, promoting the 
approach of nuclei. As noted in [4], during nuclear fusion one should expect that 
the PV is in an excited state and has an energy density comparable to nuclear 
energy. This means that the frequency and energy of PV oscillations, which 
include the appearance of (e-e+) pairs, are high. 


It is possible that when nuclei approach each other at distances d < 10 “ m, 
according to a similar scenario, an electromagnetic contribution to the attraction of 
nuclei occurs due to (x * =”) - pairs when the x” - meson is localized in the central 
region of the CC activator. 


It is interesting to note that an interpretation of the complete reduction of 
short-range strong interaction to electromagnetic interaction is also proposed [7]. 


It is obvious that if there is a sufficient number of (ee) - pairs in the 
composition of the CC activator, the process of bringing nuclei closer to distances 
2d, significantly smaller than R%, becomes real. Therefore, there is no need to focus 
on achieving Rj. as a prerequisite for starting the second stage of nuclear fusion. 
This, in turn, ensures the start of strong nuclear interaction, removing the question 
of overcoming the Coulomb barrier. 
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Abstract 

Traditional thermodynamics explains what is witnessed when 
heating matter, at both low and high temperatures, in terms of 
changes in heat capacities. To enhance our understanding, what 
constitutes thermal energy will be analyzed. It will then be 
demonstrated that it is actually the relationship between thermal 
energy density and temperature that changes, thus allowing heat 
capacities to theoretically remain constants. 
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1. Introduction 


For: 400 K > T > 250 K, condensed matter has a heat capacity (C,) that 
can be approximated by: 


Cy-souia © Cy-tiquia © 3R = 24.93 J/(K mole) (1) 


Where the subscript “v” signifies the process being isometric (A.K.A. 
isochoric: dV=0) and “R” is the ideal gas constant: R = 8.31 J/mole. 


Although not as clear in Fig. 1, for: 400 K > T > 250 K, a monatomic 
gas has an isometric heat capacity that can be approximated by: 


Cy-gas ¥ 3R/2 = 12.46 JK mole) (2) 


Similarly, for: 400 K > T > 250 K, a polyatomic gas has its isometric 
heat capacity approximated by: 


Cy-gas = R(n" + 1/2) (3) 


where “n’* represents the number of atoms in a gaseous molecule. Note 
that eq. 2 and 3 have been explained in terms of this author’s improved 


kinetic theory [2]-[4]. 


When an isobaric system expands, it performs lost work onto/into the 
surrounding atmosphere, as defined by: 


WMrost = (PAV) atm = nRdT (4) 
Lost work [(PdV) atm] has been described by this author as the energy that 
is expended by an expanding system when lifting the overlying 


atmosphere’s mass. Ultimately, it signifies an atmospheric potential energy 
increase [4]-[9]. 


Therefore, for: 400 K > T > 250 K , an expanding monatomic gas has an 
isobaric (dP=0) heat capacity that can be approximated by [4],[5]: 


Cy-gas = Cy-gas + R = 5R/2 = 20.78 JK mole) (5) 


where the subscript “p” signifies the process as being isobaric (dP=0). 
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Similarly, for: 400 K > T > 250 K, a polyatomic gas has an isobaric heat 
capacity that can be approximated by [2]-[4]: 


Cy-gas = R(n" + 1/2) +R = R(n" + 3/2) (6) 


Equations 1, 2,3, 5 and 6 are all valid approximations. They are focused 
upon the molecule’s motions. They do not fully consider all the EM 
components in either gases or condensed matter. Furthermore, equations 
2,3, 5 and 6 do not take into account any energy associated with any photons 
residing in the freespace surrounding the gas molecules, 


Furthermore, eq. | is isometric. If the condensed matter expands in some 
isobaric process then one should consider any work done in lifting the 
atmosphere, as defined by: 


Wrst = (P AV atm (7) 


As illustrated in Fig. 1, a solid’s heat capacity asymptotically approaches 
zero, as the solid’s temperature approaches absolute zero (T > 0). 


The accepted theory for heat capacities when approaching absolute zero 
has a certain basis in statistical thermodynamics. A theory that originates in 
the late 19"/early 20" century from the likes of Maxwell and Boltzmann. 
This theory is also partially founded in Einstein’s (1907) photoelectron 
model which is based upon a 3-d crystalline solid possessing a large 
number of quantum harmonic oscillators. Einstein’s theory qualitatively 
explained the heat capacities as T — 0, but quantitatively there are problems 
with the analysis. This seemingly was resolved with Debye’s analysis. ! 


Debye (1912) extended Einstein’s theory with the solid-state equivalent 
of Planck’s law for blackbody radiation, which treated EM radiation as a 


’ As a side note, Einstein’s theory can be taken as one of the early quantum-founded 
analysis. Interesting, Einstein, Podosky and Rosen (EPR argument) eventually argued that 
quantum theory is an incomplete theory [10]. Although quantum theory is currently 
embraced by the majority, some still argue that there are unresolved issues. One issue being 
that energy releasing processes tend to be irreversible, while quantum is limited to 
reversible phenomena [11]. Einstein further noted that quantum mechanics can only 
represent point-like particles in a vacuum. Of interest, this author’s understandings of lost 
work, limits reversibility of expanding systems to vacuums [4], something that entropy- 
based traditional thermodynamics fails to elucidate, i.e., work must be done onto something 
e.g., lifting of the overlying atmosphere’s mass against Earth’s gravity by an expanding 
system [4]-[9]. 
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photon gas confined in a vacuum space. Debye 7° law arrived is accepted 
as being a superior fit to experimental findings. 


Debye’s 7° law considers that the isometric heat capacity at low 
temperatures, varies as the cube of temperature, e.g.,: 


Cy-souia = A"(T/0)*) = AT? (8) 
where “A” ” and “A” are constants and 8 is the Debye temperature. 


Furthermore, for an ideal gas, the accepted theory for molar heat 
capacities of gases is based upon the degrees of freedom (/) analysis, i.e.: 


c= (Oe ° 


and 


Cp =Cy+R=(E)R+R (10) 


Degrees of freedom arguments are founded upon equipartition where a 
gas molecule has a mean energy of: kT7/2 per degree of freedom, with “k” 
being Boltzmann’s constant. Such degree of freedom arguments should be 
questioned because they assert that a gas molecule automatically obtains a 
mean energy (KT/2) per degree, without any clarity explaining the energy’s 
origins. 


Furthermore, traditionally one relies upon quantum arguments, e.g., 
cooling molecules move from a higher to a lower quantized energy state. 
Moreover, the fraction of molecules in the lowest energy state approaches 
unity, as the temperature approaches absolute zero. Such arguments have 
their place when contemplating the excitations of electrons but thermal 
energy may be better contemplated in terms of broad continuous spectrums 
of energy rather than discrete/quantized energy states. 


A goal of this paper is to explain what has been witnessed without relying 
upon Einstein’s analysis, or Debye’s 7° Law. A result will be an improved 
understanding concerning what constitutes thermal energy. A result that 
does not rely upon discrete quantized states. 


2. The Essentials For: 400 K > T> 250 K 


Eq. 1 can be employed to approximate the thermal energy (Qin) required 
to raise the temperature of m moles of condensed matter by dT: 


Qin-v = 3NRAT = NCy_soyqgdT = NCy-tiquiadT (11) 
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For the isometric heating of » moles of a monatomic gas, eq. 2 can be 
used to approximate the thermal energy (Qin) required to raise the gas’ 
temperature by dT: 


Qin-v = (3) nRAT = nCy-gasdT (12) 


Similarly, for the isometric heating of m moles of a polyatomic gas, eq. 3 
can be used to approximate the thermal energy (Qin) required to raise the 
gas’ temperature by dT: 


Qin-v = (n" +1/2)nRdT = nCy_gasdT (13) 


And, for the isobaric heating of m moles of an expanding monatomic gas, 
eq. 5 can be used to approximate the thermal energy (Qin) required to raise 
the gas’ temperature by dT: 


Qin-v = (n" + 3/2)nRAT = NCy-gasdT (14) 


Similarly, for the isobaric heating of n moles of a polyatomic gas eq. 6 
can be used to approximate the thermal energy (Qin) required to raise the 
gas’ temperature by dT: 


Qin-p = (2) (n" + 3/2)mRAT = NCp_gasdT (15) 


For: T < 250 K. Conversely, for the heating of n moles of condensed 
matter, eq. 8 can be used to approximate the thermal energy (Qin) required 
to raise the temperature by dT: 


Qin-v =ATFAT = Cy-souadT (16) 


Equations 1 through 6, eq. 8, and equations 11 through 16, all assume 
that the heat capacities are constants and that the thermal energy is some 
linear increasing function of temperature, i.e., for T> 250 K: 


aT = T; — T; (17) 
where the subscripts “/’ and “i” 
states. 


3. The Essentials For: T< 250K 


For: T < 250 K, one can integrate eq. 16 (based upon Debye’s T? law), 
thus obtaining for condensed matter: 


Qin-y = Sddin-y =AS,'T3dT = (A/4\(TH-T) (18) 


i 


respectively signify the final and initial 
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Generally condensed matter can be contemplated as being isometric. 
Therefore, one can contemplate eq. 18 in terms of energy density, thus: 


Qin-v/V = (A/4V) (Tf -T/) (19) 


Are there any relations in thermodynamics that are a function of 7%? 
What comes to mind is the Stefan-Boltzmann’s equation for a blackbody 
radiation spectrum: 


pp = oT* (20) 
where “pz” is the blackbody radiance (energy density), and “oa” is the 
blackbody radiation constant: 0 = 872°k*/15c?h? 


= 5.670374419...x10°° (W-m?-K), 

The constants in eq. 19 and eq. 20 are not one and the same, ie., a # 
(5). However, the blackbody energy flux as described by eq. 20 is founded 
upon the thermal energy emitted by certain closed system’s walls, into that 
closed system’s freespace. 

4. The Essentials, For All 7 


Certainly, all condensed matter above absolute zero emits energy whose 
broad continuous spectrum often approximates a blackbody spectrum. This 
radiated thermal energy is a part of what can be referred to as radiant energy. 
Radiant energy can take many forms ranging from gamma rays to 
microwaves. The question then becomes, what part of the spectrum is 
related to thermal energy in condensed matter? 


In order to answer that question, one starts by considering the energy 
density [u(v, T)] associated with blackbody radiation in a closed system. As 
defined by Planck’s formula: 


w(v,T) = (22) vr jeter (21) 


The wavelength version of eq. 21, for the energy density of a blackbody 
spectrum between A and d/ is the energy per mode: (EF = “ = hv), times 


the density of states per photon, times the probability that the mode is 
occupied. That being: 


u(A)da = (=E) (esa (22) 


Plank’s eq. 22 is useful for calculating the total energy density (J/m*) 
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associated with a full blackbody spectrum in an enclosed freespace. Based 
upon eq. 22, plotting the total energy density, one obtains Fig. 2. 


Let’s reconsider condensed matter. From an energy density perspective 
condensed matter concentrates thermal energy, when compared to the 
surrounding radiation density in freespace. It is inarguable that when 
considering a vacuum in thermal equilibrium that the vacuum’s walls must 
absorb as much thermal energy as they emit. 


The exact mechanism for the emittance and absorption of thermal 
spectrum lacks any clarity. Interestingly, this author has previously 
proposed that such thermal spectrum is due to inelastic collisions (both 
intermolecular and intramolecular), Furthermore, a plausible explanation 
concerning the mechanism that creates the thermal photons associated with 
radiated thermal energy may be the physical distortions of the electron 
clouds within molecules and/or atoms, during inelastic collisions [4], [12]. 


Thermal photons are considered to belong to a broad continuous EM 
spectrum. The part of the EM spectrum of energy that tends to be both 
absorbed and emitted by matter. 


Whatever the mechanism turns out to be, it must explain why thermal 
spectrums are: 


1) Universal to all condensed matter above absolute zero. 
2) Have a temperature functionality. 
‘3) Disappears as: T > 0 K 


Interestingly, the above 1), 2) and 3) also applies to the thermal energy 
density within matter. However, the temperature functionality for matter 
energy density is approximately linear, at least for: 450 K > 7 > 250 K. 
Conversely, a blackbody’s energy density clearly is not linear, as is 
illustrated in Fig. 2 


Strangely, for T < 250 K, based upon Deby’s energy density (eq. 19), 
one would get a graph similar to that depicted in Fig. 2 


As previously stated the above has traditionally been reconciled by 
contemplating that heat capacities are temperature dependent, rather than 
constant through all temperature regimes. 


5. New Graphical Analysis 


The problem with eq. 22 (hence, Fig. 2) is that it considers all 
wavelengths of a blackbody spectrum. Thermal energy only concerns the 
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wavelengths (frequencies) that tend to be absorbed by matter, thus 
becoming thermal energy within that matter. It is accepted (without absolute 
clarity) that wavelengths associated with thermal energy predominately 
resides in the infrared part of the spectrum. Therefore, contemplating the 
total energy density (full spectrum), as sketched in Fig. 2, can be 
disingenuous. 


Accepting that thermal energy consists of only part of the total energy 
density. Then in order to calculate the thermal energy density, one only 
needs to consider the part of the spectrum that is predominately thermal 
energy. What parts of the EM spectrum does one quantify as being thermal 
energy? 


The thermal spectrum resides in the infrared: 0.75 um = 750nm < A< 
100 wm = 100,000 nm However, the division between visible light and 
infrared is founded upon our human vision rather than any concise notions 
of thermal energy. Accordingly, this analysis depends upon where one 
chooses to start their consideration of what constitutes thermal energy. 


In remote sensing (often used by geologists) the thermal spectrum is 
considered to be in the range of: 3 um < A < 15 wm. In thermodynamics 
the thermal infrared is considered to be in the range of: 8um<A< 
15 um. The difference can be attributed to the fact that remote sensing 
sometimes involves spotting forest fires whose emitted wavelengths are 
close to 3 wm. Again, there is no universal understanding concerning 
thermal energy. 


Based upon Fig. 1, one can speculate that the transition between linear 
and T? functionality for the energy density within condensed matter, occurs 
near 250 K. Thus, one begins the analysis by examining the energy density 
for: T = 250 K, as illustrated in Fig. 3. One immediately notices that the 
peak is near 11.4 wm. Consider that the thermal energy starts somewhere in 
the range: 8um <A<12 um. 


Unfortunately, the integration (of eq. 21 and/or eq. 22) can only be 
performed over all wavelengths (0 < A < &) (this is accomplished using 
mathematical tricks). Therefore, one cannot simply calculate the thermal 
energy density associated with only some part of a blackbody spectrum. 
However, one can graphically approximate the thermal energy density 
associated with various parts of the blackbody spectrum. 
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One starts by sketching and graphically calculating the total energy 
density for all temperatures between 50 K and 600 K degrees. This was done 
at 25 K intervals using graphs obtained from a website that provides the 
“Planck Curve for Blackbody Radiation”, as was published by TU Graz 
University of Technology [13]. 


The total energy density of blackbody radiation was obtained using 
geometric approximations. For example, adding the areas of triangles A, B 
and C plus the rectangle D, as illustrated in Fig 4 for four different 
temperatures. The result for the total thermal energy was plotted in Fig. 2. 


A similar graphical technique was then used to calculate what is then to 
be considered as the “thermal energy density”. A cut-off wavelength (A,,,¢) 
between what is considered thermal energy and non-thermal energy was 
determined. 


Three cut-off wavelengths were chosen: Apy,; = 8 um, Acyz = 10 um 
and A, = 12 um. Fig. 5, illustrates the thermal energy density at four 
different temperatures (100, 200, 300 & 400 K), for Acyy = 10 um, hence 
over the range: 10um < A < ©, 


The data table for total energy density and thermal energy density, as 
well as the energy flux (based upon the Stefan-Boltzmann equation: eq. 20) 
is given in Table 1, at the end of the paper. 


Plotting the total energy density, and the total thermal energy density 
versus temperature for the three cut-off wavelengths (A,,,;), one obtains Fig. 
6. It is of interest that the total thermal energy density can be approximated 
as a linear relation for: 250 K < T< 400 K. The approximate slopes for the 
thermal energy density using the three cut-off wavelengths are: 


Deut = 8 Um: (12.58-2.59/ (400-250) = 6.6 x 107 m/K 
Acut = 10 pm: (9.45 -2.27)/(400-250) = 4.8 x 10? J/m3/K 
eur = 12 um : (6.99 — 1.95)/(400-250) = 3.4 x 10? J/m3/K 


It must be emphasized that this analysis was formulated to show that, if 
the thermal energy density had a thermal energy cut-off wavelength, then 
one may be able to explain why thermal relationships are approximated as 
being linear functions of temperature for: 250 K < T< 400K. 


Furthermore, for: T < 250 K the thermal energy density would be a 
rapidly decreasing function of temperature. Thus, it becomes likely that the 
heat capacity approximates being a constant through all temperatures. 
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Moreover, it is the thermal density that is not a linear function of 
temperature over all temperature regimes. 


It becomes a case of our residing on Earth’s surface, that we have 
intuitively come to falsely believe that thermal energy is a linear function 
of temperature. 


In accepting this analysis, one assumes that the thermal energy density 
associated with matter’s thermal spectrum, can be related to the thermal 
energy density within condensed matter. Based upon the fact that they are 
both linear relations of temperature, one must conclude that such a 
relationship may exist. Unfortunately, at this point this relationship is not 
fully realized. Of interest, one should include the fact that the thermal 
spectrums emitted by condensed matter, tend to resemble cut-off blackbody 
spectrums. 


Furthermore, the thermal energy cut-off at: 8 um, 10 um, and 12 wm, 
were somewhat arbitrarily made. Certainly, a superior cutoff may be 
obtained by performing a more exacting analysis over a broader range of 
wavelengths and by considering various types of condensed matter. Note 
that this author also tried a cutoff of 3 um with a preliminary determination 
that it was too low (see the next section: Weakness of Graphical Analysis). 


Moreover, the cut-offs were chosen to be vertical and linear, when in 
fact it may vary and actually be dependent upon the nature of the condensed 
matter itself. Even so, there should be a certain degree of universality. 


A point becomes that one does not necessarily have to rely upon 
Einstein’s quantum arguments and/or Debye’s 7° law to explain why the 
total energy required to heat matter has temperature dependence. Again, the 
simpler explanation may ultimately reside with the realization that it is the 
thermal energy density that changes with temperature rather than the heat 
capacity varying in different temperature regimes. 


6. A Weakness in this Graphical Analysis 


Looking at Fig. 5 one will notice that as the temperature increases, the 
cut-off wavelength of interest shifts right with increasing temperature. The 
reason that this analysis was stopped at 7 = 500 K is that the cut-off 
wavelength approached C in Fig. 5. At which point, this graphical analysis 
becomes increasingly inaccurate (see the data points with ** in Table 1). 


The expectation should be that as the temperature increases, the thermal 
energy density continually increases, which is not necessarily what one gets 
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when doing simple time- consuming geometric analysis, as was done 
herein. It would be of interest for someone to extrapolate this analysis, 
perhaps using computing, e.g., artificial intelligence (AI). 


7. What does the Cut-Off Represent? 


Whether it is ultimately placed at 8 wm, 10 wm, 12 wm, or some other 
value, the cut-off wavelength (A,,) represents a theoretical limit for the 
thermal energy, that being the energy associated with the kinematics within 
matter. A result of such kinematics is energy radiating from all condensed 
matter that is above absolute zero. With kinematic energy being the kinetic 
(translational plus rotational) plus vibrational energy in matter. This 
statement is partially based upon this author’s kinetic theory [2]-[4] wherein 
the separation of these energies differs from the traditionally accepted 
degrees of freedom, based theorization. 


As stated, there are various types of radiant energy. These include 
microwave, visible light, ultraviolet, x-rays and gamma rays. What needs to 
be added to this list is thermal energy, which should be placed between the 
microwave and visible light spectrums. Herein, the proposal being that 
thermal energy now starts at the cut-off wavelength and extends into longer 
wavelengths, e.g., the far infrared, and perhaps even into the microwave. 


Again, for condensed matter, thermal energy includes the energy 
spectrum associated with the vibrational energies inside that matter. It 
represents the wavelengths/frequencies of energy that are radiated and 
(most likely) absorbed by that matter. 


Importantly, thermal energy is not limited to the discrete energies that 
are readily absorbed and radiated by electrons, i.e., the quantized energies 
often determined in infrared spectroscopy. When contemplating closed 
systems, thermal energy is the energy that is both radiated and absorbed by 
the walls, which often can be blackbody in nature. 


Thermal energy is a broad continuous spectrum of energy. Is it plausible 
that such a broad spectrum is somehow absorbed by the atoms and/or 
molecules, as a whole? That would certainly help one to understand why 
thermal equations involve both the number of molecules and/or the number 
of atoms. 


A possible implication being that this energy is absorbed by nuclei and/or 
their components. This plausibility has previously been briefly pointed out 
by this author [4], [12]. Whatever the exact mechanism for the absorption 
and emission of thermal energy, its result is the vibrational energy within 
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condensed matter. Importantly, since the heat capacities are based upon the 
number of atoms in molecules, there must be some sort of approximate 
universality. A universality not readily explained by either the configuration 
of its electrons, or the shapes of their electron clouds. 


8. Gases 
Can any of the above be applied to gases? 


For monatomic gases the heat capacities have been deemed to be 
relatively constant. However, that is not the case for polyatomic gases over 
all temperature regimes. 


Fig. 7 shows the isometric specific heat of air (99%: Oz and N2) over a 
broad temperature range: 0 < T < 2000 K. Note that engineers prefer to use 
specific heat which is related to the gas’s mass. Physicists and chemists 
prefer to use molar heat capacities, which are related to the number of moles 
of gas. 


As shown in Fig. 7, for, T< 200 K, a polyatomic gas’ isometric specific 
heat [heat capacity per unit mass (cy)] has been deemed to be some rapidly 
increasing functions of decreasing temperature, e.g., as: T > 0 K. 


Furthermore for, 200 K < T<375 K, apolyatomic gas’ isometric specific 
heat (cy) has been considered to approximately be constant, as described by 
eq. 3. 


Furthermore, for, 375 K < T < 1750 K, a polyatomic gas’ isometric 
specific (cy) has been considered as being a linear increasing function of 
temperature. 


Note that Fig. 7 extends to 2000 K, while the analysis for solids stopped 
at 500 K. This was due to the previously stated weakness in this graphical 
analysis. 


Interestingly, the fact that monatomic gases have a relatively constant 
heat capacity as defined by eq. 2, implies that not all gases behave the same. 
Seemingly, gases that have vibrational energy as part of their total thermal 
energy are the only gas’s whose heat capacities have temperature 
dependence. That being the case, then the kinetic energy (translational plus 
rotational) [2]-[4] is not affected, at least in the same way that vibrational 
energy is. 


Is there anything that confirms the above? An indirect proof resides with 
the density of gases being a function of temperature. Specifically, the mean 
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molecular volume of a gas should be related to its mean kinetic energy: 
3kT/2. Ditto, a gas’ mean molar volume is related to its mean molar kinetic 
energy: 3RT7/2. In other words, as the gas’ temperature increases its kinetic 
energy increases, Therefore, the gas’ mean molecular volume should 
increase, 


Conversely, as the gas’ temperature decreases the number of moles of a 
gas per unit volume (m*) should increase. As shown in Fig. 8, the number 
of moles of a gas per unit volume, increases linearly as the temperature 
decreases from 320 K to 160 K. As the temperature decreases below 160 K, 
i.e., T > 0 K, the number of moles per cubic meter increases rapidly. 


The above is expected. That is if the gas’ thermal energy density 
decreases rapidly as: T — 0 K. Perhaps, what was not expected is this. As 
shown in Fig. 8, this equally applies to both monatomic and diatomic gases, 
namely hydrogen (H2), helium (He), nitrogen (N2), argon (Ar) and oxygen 
(O2). The implication being that both monatomic and diatomic gases 
witness the same degradation of kinetic energy, as: T > 0 K. 


An explanation for the similar degradation of kinetic energy is as 
follows. This author’s kinetic theory implies that the larger structured wall 
molecules impose their kinematics upon the smaller gas molecules [2]-[4]. 
This means that both the monatomic and diatomic gases will obtain the 
same mean kinetic energy from their collisions with the structured vibrating 
wall molecules, as is defined by: 3k7/2. 


The overall implication is that it is the changes to a polyatomic gas’s 
absorption of vibrational energy that separates the energetics of polyatomic 
gases from monatomic gases, as: T > 0 K. 


9. Thermal Energy Density of High T Gases 


For a diatomic gas (e.g., 99% of air), eq. 3.can be used to approximate 
the isometric heat capacity: Cy_gqs * 5R/2, which it clearly is for: T= 250 
K (see, Fig. 7) and remains approximately so for: 250 K < T < 375 K. 
Unexpected is that its heat capacity then increases with temperature. It 
approaches the theoretical isobaric heat capacity for a diatomic gas, as 
defined by eq. 8. That being C,_ga, ~ 7R/2 at high temperatures (blast 
furnace temperatures, e.g.: 1,500 K) 


How does one explain the seemingly heat capacity’s increase for: 375 
K <T<1750 K, shown in Fig. 7? Consider the heat capacity to be constant, 
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then again, it is the required thermal energy in per degree temperature 
change that increases. In other words, Qin-v as defined in eq. 13 increases 
because there is an unexpected increase to the thermal energy density in 
relation to the gas’s temperature increase. Obviously, the gas’ thermal 
energy density increase is greater than eq. 13 implies, as illustrated in Fig. 
7 for: T > 375 K. 


An implication of Fig. 7 is that the thermal energy density of the 
sutrounding thermal (or blackbody) radiation increases rather rapidly 
between: 450 K< T < 750 K. After which, the increase tapers off: 750 K < 
T < 1750 K. Seemingly approaching some linear increasing value as 
temperatures exceed 2000 K. It would be of interest to attain a superior 
analysis of the thermal energy density at higher temperatures than has been 
provided by the graphical analysis herein. Again, hopefully such an analysis 
may be improved, e.g., by using AI computing. 


10. Thermal Energy Density of Low T Gases 


The traditional explanation for the heat capacity of a gas increasing as 
temperatures approach absolute zero, (i.e. for: T< 250 K, as: T — 0) was 
based upon the concept that states can be frozen-out. Specifically, the 
vibrational energy of gases (e.g., air), are increasingly frozen-out as the 
temperatures decrease. Perhaps, but air is 99% homonuclear O2 and N2, both 
of which traditionally speaking, are supposedly transparent to thermal 
radiation. 


Imagine that O2 and Nd truly are transparent to thermal 
radiation/photons, as has been traditionally claimed. One, might then ask, 
what difference would being frozen-out make? One faces that enigma 
bestowed in 20" century thermodynamics. That being, energy does not 
simply exist, as is implied by equipartition and degrees of freedom, founded 
kinetic theory. 


Energy associated with a gas molecule must have an origin. For the case 
of a diatomic gas molecule in a closed system, this author’s kinetic theory’s 
origins are [2]-[4]: 


1) Kinetic energy (translational plus rotational) is primarily imposed 
onto the diatomic molecule through the gas molecule’s collisions 
with structured vibrating wall molecules. 

2) Vibrational energy is primary obtained from the absorption of 
thermal photons. 


Additionally, if homonuclear diatomic gases (e.g., O2 and N2) actually 
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are transparent to all thermal radiation, then being frozen-out would have 
no consequences. In which case, the diatomic gas must obtain both its 
kinetic and vibrational energy from its collisions with walls. Based upon 
heat capacities, this seemingly is not the case [2]-[4]. 


Importantly, thermal energy consists of broad continuous spectrums of 
energy, e.g., blackbody spectrums. Again, this differentiates thermal energy 
from quantized energies associated with the discrete energy states that are 
associated with electrons. It raises the question, is there another plausible 
explanation? 


The simplest explanation for what is witnessed for gases at low 
temperatures is that the thermal energy density rapidly approaches zero, as: 
T > 0K. Thus, one has the illusion that the heat capacity has skyrocketed, 
when in fact it is thermal energy density that has crashed. 


In terms of polyatomic gases (e.g., air: 99%: N2 and O2) any thermal 
energy absorption resulting in vibrational energy would rapidly approach 
zero as: T ~ 0 K. Hence, the vibrational energy of that gas would also 
rapidly approach zero. 


One may prefer to think in terms of the gas’ vibrational energy being 
frozen-out at low temperatures. In many ways both arguments have similar 
results. A point becomes that one is not restricted to quantum arguments in 
order to comprehend why air’s heat capacity has the illusion of being some 
rapidly increasing function of temperature as: T > 0. 


Moreover, thermal energy concerns the motions and vibrations of 
atoms/molecules as wholes. The energy associated with such motions 
quickly becomes zero as: T ~ 0K. This is not the same as the energy 
associated with electrons. It is one thing to say that molecular motions 
cease, however this does not imply that electrons necessarily stop moving. 


11. O2 and N2 Are Not Completely Transparent 


How can diatomic gases like O2 and N2 be transparent to thermal 
radiation, yet their heat capacities adhere to standardized equations 2, 3, 5 
and/or 6? 


Heat capacity experiments involve the measurements of thermal energy 
in (Qin ) versus the temperature change (dT). The required energy is based 
upon the gas’s kinematics, i.e., the gas’s kinetic plus vibrational energies. 
For all intents and purposes, thermal energy constitutes wavelengths that 
are greater than the cut-off wavelength. That being the previously described 
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theoretical division for thermal energy, i.e., the thermal part of the infrared 
spectrum. 


The concept of certain gases being transparent is founded upon infrared 
spectroscopy, which often focuses upon shorter wavelengths photons 
(higher frequencies/energies). Those, being the photons that tend to interact 
(as waves) with various gases electron configurations, e.g., a gas’ lopsided 
charge distribution. To what extent the absorption due to different electron 
configurations actually results in the heating of a gas, may become 
debatable on a case-by-case basis. 


Many believe that the absorption of discrete wavelengths due to electron 
configurations is the primary mechanism by which gas molecules absorb 
thermal energy. In terms of this paper, so-called transparent gases likely 
absorb thermal photons. They do not absorb the higher energy (shorter 
wavelength) photons that tend to excite discrete energy states in gases like 
CO2. Note that CO2 also has a strong absorption peak near 15 um. 


It is of interest that in remote sensing, the wavelengths in the thermal 
part of the EM spectrum forms one of the windows for observing the Earth 
from satellites. Seemingly, the majority (up to 80%) of the photons in this 
window do readily pass through the atmosphere. An implication being that 
the absorption of thermal photons by gases is inconsequential. 


Due to the speed of light, the absorption of thermal photons by gases can 
represent relevant quantities of thermal energy. Imagine a closed system 
where thermal photons are both radiated and absorbed by the walls. 
Certainly, the chances of a given photon scattering with and then being 
absorbed by a given gas molecule, remains low. However, the quantity of 
photons that can interact with that gas molecule over a given duration 
renders the absorption of thermal photons by gas molecules, real. 


12. Infrared Spectroscopy 


Would all of a gas’ absorptions be witnessed in infrared spectroscopy? 
In dispersive infrared spectroscopy, the resulting sample’s spectrum is 
obtained by subtracting a background (or zero) spectrum from one’s 
measurements. The resultant spectrum.then shows any uniqueness which 
often is based upon the sample gas’s unique electron configuration, e.g., 
lopsided charge distributions. This uniqueness enables one to quickly 
identify the gas. 


Consider that Fig. 9. Path 1 (reference cell path) will have the same 
interference (absorptions and radiances), as Path 2 (sample cell path). The 
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interference along Path 1 when measured gives the background spectrum. 
This interference comes from the air, contents of the reference cell, 
laboratory and/or apparatus etc. The subtraction of the background 
spectrum means that the recorded spectrum only shows uniqueness 
associated with the sample being tested. A uniqueness that is based upon 
discrete energy associated with specific wavelengths, which are amplified 
before the spectrum is recorded. 


It should be emphasized that both the analyte cell (sample cell) and the 
reference cell are enclosures. The reference cell is either evacuated 
(vacuum) or is filled with a gas that is considered as being non-absorbing, 
e.g., air (99% homonuclear O2 and N2). As enclosures, both cells contain 
thermal (e.g., blackbody) radiation that is related to their temperature. 
Furthermore, the light source is a heated wire radiating a thermal (e.g., 
blackbody) spectrum that is associated with its higher temperature. 
Therefore, much of the subtracted background will constitute a broad 
continuous thermal (e.g., blackbody) spectrum, whose peak is in the vicinity 
of the previously described cut-off wavelengths, e.g., Agu: = 10 um. 


Few have asked, does the background subtraction have any 
consequences, other than enabling the identification of the gas? If the 
background spectrum includes some universal phenomenon, then the 
resulting spectrum becomes misleading from a thermal energy perspective. 
For example, if both the analyte and the reference gas absorbed (or emitted) 
the similar parts of the light source’s spectrum. Then it becomes erroneous 
to assume that the parts of the subtracted background spectrum are in no 
way universal to all (or most all) gases. In which case, the background 
subtraction includes the subtraction of certain gas’ thermal signature, or part 
thereof. 


Similar argument holds for Fourier Transform Infrared Spectrometry 
(FTIR), where the background’s subtraction occurs in its computer and 
programming. Note that in FTIR the background spectrum is often 
measured in the cell prior to testing the sample. Its data is stored in the 
computer, where the background spectrum is subtracted. 


Interestingly, in inelastic collision studies, infrared spectroscopy is used. 
Furthermore, as previously stated most such studies are performed in 
ultracold experimental apparatus, hence certain aspects of the thermal 
spectrum would be eradicated. 


The plausibility that some universal phenomena may not be witnessed in 
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the infrared spectroscopy, raises concerns that may help one’s 
understanding. 


13. Discussion 


Nothing was proven here beyond all doubt. Ditto can be said for 
Einstein’s photoelectron model which treats solids as a group of quantum 
oscillators. The same applies to Debye’s 7° law. Certainly, both Einstein’s 
and Debye’s models are based upon mathematical results. Results that can 
be tied to both statistical and quantum theory. 


The issue with such results, is that one often places mathematics ahead 
of theoretical logic, i.e., one’s pet theory is purely mathematical thus, 
leaving theoretical justification by the wayside. In this paper a theory was 
devised ahead of the math, a theory whose plausibility was verified using 
graphical analysis. A theory that fits with and clearly explains accepted 
thermodynamic relations. 


It is interesting to note that both Stefan-Boltzmann’s and Planck’s energy 
density relations, are in many ways part of quantum theory’s origins. 
Interestingly, this is the point that this theory separates from the traditional 
understandings. It raises questions, e.g., are such accepted relations founded 
upon mathematical tricks, or divine logic? Planck thought that the 
quantization of energy was a mathematical trick enabling one to advert the 
ultraviolet catastrophe [16]. Certainly, at that point Planck was thinking in 
terms of classical oscillators. 


The issue is not with the quantization of a photon’s energy, rather it is 
the notion that thermal energy is only adsorbed by electrons possessing 
discrete quantized energy states. The states that are found in infrared 
spectroscopy results, enabling one to identify specific gases. Energy states 
that when absorbed may result in thermal energy increases but equally may 
result in the emission of a photon of that discrete energy. 


Thermal energy has a universality based upon the number of atoms in a 
molecule rather than some configuration of electrons. Again, this 
universality is witnessed in equations 1 through 6 and 11 through 15. 
Furthermore, thermal energy concerns broad continuous spectrums, 
spectrums whose total energy density can be best inferred by heat capacity 
experiments. 


In terms of eq. 22, the energy density of a blackbody spectrum between 
A and dA is the energy per mode, times the density of states per photon, 
times the probability that the mode is occupied. It remains valid as a 
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mathematical approximation. It is just that the “density of states” must 
rapidly decrease at the cut-off wavelength for shorter wavelengths (for 
higher energy photons/phonons), when contemplating thermal energy 
densities. 


Again, the cut-off wavelength is likely not a sharp vertical divide as was 
graphically depicted herein. It must be acknowledged that the cut-off 
frequency may vary with states (solids, liquids and gases) as well as with 
the nature of the matter. Interestingly, some matter’s absorption of thermal 
energy includes visible wavelengths, and the near infrared. These 
absorptions are primarily due to electrons absorbing energy (as waves), 
rather than the molecules/atoms themselves. Furthermore, for gases there 
should be a pressure dependence, which affects the mean free path of the 
thermal photons. 


It is of further interest that one of our estimated cut-off wavelengths 
(Acut = 10 um) is at the end of what is referred to as the “fingerprint zone” 
(2 um < A < 10 um). That being where there are many complex bands that 
form a molecule’s fingerprint. The chosen cut-off also resides in what is 
often considered as being in the thermal infrared (8 um < A < 15 um), 


This author previously theorized that thermal radiation, e.g., blackbody 
radiation, may be due to inelastic collisions between atoms and/or 
molecules. Moreover, it was speculated that the emission of thermal 
photons was due to collision induced distortions of the electron clouds in 
both atoms and molecules [4], [12]. Thus, explaining the broad continuous 
nature of thermal spectrums, at least in the thermal infrared. 


Logic implies that if matter radiates thermal photons then that matter 
should absorb thermal photons. Questions then arise concerning what is the 
mechanism of absorption? Is it by the electron clouds, or are thermal 
photons absorbed by another mechanism. 


This author has previously discussed that absorption of thermal photons 
may involve photons behaving as particles thus passing their momentum 
onto the atoms/molecules [4], [12]. Thus, there may be more to the wave- 
particle duality of photons than has been accepted. 


Whatever the mechanism, thermal equilibrium implies that matter 
absorbs as much thermal energy as it relinquishes (kinematic plus radiated 
thermal energy). This is certainly witnessed in enclosures where the walls 
must absorb as well as radiate thermal photons. 


Accordingly, the mechanism of absorption and emission do not 
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necessarily have to be identical. The only criterion is that when in thermal 
equilibrium, the total thermal energy in equals the total thermal energy out. 
Although speculative, one could envision matter absorbing thermal photons 
as particles and emitting thermal photons due to inelastic collisions, e.g., 
collision induced distortions of electron clouds. 


It should be stated that statistical thermodynamics describes all of a 
system’s energy solely in terms of that system’s kinematic energies. Such a 
description remains a valid approximation for closed systems where the 
system’s total kinematic energy is significantly greater than its total 
radiative energy [4], [17]. And where the closed system’s wall both radiate 
and absorb thermal photons, thus giving the illusion of elastic collisions [4], 
[12], [17]-[19]. 


14. Conclusions 


We instinctively think in terms of temperature being a measure of 
warmth, realizing that as the temperature increases, then so too does the 
thermal energy. Our first experiences involve thermal energy being a linear 
function of temperature, as is measured in our daily lives using a 
thermometer. Therefore, when performing thermal energy experiments, one 
instinctively conforms to it being the matter’s heat capacity that varies with 
temperature. As a powerful mathematical construct, _ statistical 
thermodynamics fortified our confidence in what may be erroneous logic. 


Contemplating that it is the thermal energy density that varies with 
temperature, alters our thermodynamic perspective. A perspective that 
allows heat capacity to be considered as a constant over most temperature 
regimes. A new perspective that does explain our experimental findings in 
a simpler context. 


Specifically, experiments based upon the thermal energy into (Qin) a 
system may be better explained in terms of the system’s relatively constant 
heat capacity multiplied by its temperature change (dT). However, the 
amount of thermal energy required to incrementally increase the 
temperature (dT) is a not always a linear function of temperature, over all 
temperature regimes. 


There remains a relatively small temperature regime over which the 
thermal energy required to increase a system’s temperature remains a linear 
function of temperature. That happens to be the temperature regime in 
which we tend to reside here on Earth’s surface, i.e., 250 K < 7 < 400 K. 
Seemingly, this reality has skewed our perspective. 
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Graphs for Planck’s blackbody energy density relation were obtained at 
25 K intervals for: 50 K < T<500 K. The total energy density based upon 
the plots were graphically calculated. Cut-off wavelengths of 8 um, 10 um, 
and 12 um were then chosen based upon the accepted graph for the heat 
capacity of mercury, and the graph for the blackbody energy density at: T= 
250 K. 


All wavelengths longer than the chosen cut-off wavelength were 
considered to be part of the blackbody’s thermal energy. The total 
blackbody thermal energy density was then graphically calculated, e.g. for 
a cut-off of: Apyt = 10 um, the blackbody spectrum’s total thermal energy 
density was considered to include: 10um <A < ©, 


Plotting both the blackbody spectrum’s total energy density at various 
temperatures, it became obvious that as expected, the blackbody’s total 
energy density is a rapidly increasing function of temperature. 


However, the blackbody spectrum’s thermal energy density could be 
approximated as a linear function of temperature for: 250 K < T< 400 K for 
all three proposed cut-off wavelengths. Assuming that the blackbody 
spectrum’s thermal energy density must be related to the thermal energy 
density in matter, gives clarity as to why most of our accepted thermal 
energy relations are linear increasing functions of temperature. 


For temperatures below 250 K the thermal energy density decreases 
rapidly (e.g., T*) as the temperature decreases, i.e., as: T — 0 K. At absolute 
zero the thermal energy density is considered to be zero, i.e., all molecular 
motions cease. 


Conversely, for temperatures above 400 K the thermal energy increases 
more rapidly than a linear relation of temperature would indicate. 
Unfortunately, the graphical analysis became unsensible for: T> 500 K. 


It must be emphasized that this is a preliminary paper, written to show 
that there is a simpler logic that explains what has been witnessed. The 
following all warrant more consideration: 


e Cut-off wavelengths (e.g.,A¢y¢ = 10 wm) for thermal energy. 

e Nature of thermal energy density for: T > 500 K. 

e Relationship between a blackbody spectrum’s thermal energy 
density and thermal energy density in specific matter. 


Ultimately, this preliminary analysis separates new thermodynamics from 
the traditionally accepted theories based strictly upon statistical 
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thermodynamics, and quantum physics. 
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Fig. 1 is a sketch of the heat capacities versus temperature for mercury. 
Note that discontinuities exist where phase changes occur. Also note that 


mercury is unreactive, thus when in the gaseous state mercury tends to be a 
monatomic gas [1]. 
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Fig. 2 is a sketch for the total energy density versus temperature, as 
described by eq. 22. 
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Fig. 3 is a sketch for the blackbody energy density versus wavelength at: 
T=250K 
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Fig. 4 shows the total energy density versus wavelengths for 7=100, 
T=200, T=300 and T=400. The blackbody energy density is the summation 
of the areas in triangles A, B and C plus the rectangle D. 
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Fig. 5, shows the thermal energy density versus temperature for cut- 
off: A = 10 ym. For: T = 100 K it is the summation of: A, B, G& D. For: T 
= 200 K it is the summation of: A, B, CD & E. For: T = 300 K and T= 400 
K, itis the summation of B, C& D. 
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Fig. 6‘shows both the total energy density and the thermal energy density 
versus temperature associated with blackbody radiation for the three chosen 
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cut-off wavelengths. 


Specific Heat of Air versus Temperature 
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Fig. 7 is a sketch of the specific heat of air versus temperature. Note that: Ri7-RMa 


that being the dry air constant where Mz is the dry air mass and R is the universal - 
gas constant [14]. 
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Fig. 8 shows the number of moles per cubic meter of various gases as a 
function of temperature. Data is in Table 2 at the end of this paper. It is 


based upon data from the “Handbook of Chemistry and Physics” [15]. Note 


that this also implies that Avogadro’s hypothesis comes into question at low 
temperatures. 
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Fig. 9. Is a schematic for a dispersive infrared spectrometer. 
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Table 1: Gives Blackbody thermal energy at various temperatures 


AS Peak Energy | Thermal Energy Density 

Toya [u(v,7)] | Eq.20 | Density | At three chosen cut-offs 

(K) | (um) | (J/m3) 

toe ee 
sum | 10pm | =12 um 
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Be Naneersible due to weakness in graphical analysis 
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Table 2 


Density at 1 Bar Number of moles per m? at 1 Barr 
pressure (kg/m?) pressure 


rw fc 
omic 
mass 

H2 | He |Ar | No | O2 ; 


| Oo 132 [250] | | {656.25 | 624.19 | | 
| 40 [062 [1.20] [| jf 30541 [29985 | 

[041 [oso] | fj 20109 | 20005, | TT 

| 0.30 [0.60] [| 438 | {| 150.55_ | 150.10 | | 156.31 | 
| 0.40 | 4.06 | 2.84 | 3.25[ 100.20 | 100.12 | 101.58 | 101.38 | 101.63 _| 
|_280_| 0.09 | 0.17] 1.72 | 1.20 11.38{ 42.93 | 42.93 | 42.98 | 42.98 | 43.00_| 
| 300 | 0.08 [0.16] 1.60] 1.12 [128] 40.06 | 40.08 | 40.13 | 40.09 | 40.13 | 


Handbook of Chemistry and Physics: Editor: W.M. Haynes (CRC Press, 
ew York, 2010) [14] 


| O11 | 
| 0.10 _| 
[0.09 _| 


| 0.19 | 


AW signifies Atomic Weight. Note: mols/m’ calculated by: 
Density* 1000/A W 
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